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he  following  improvements  have  been  made: 

\>  f%  Composite  materials  of  laminated  filamentary  construction  can  be 
resized  for  strength  and  flutter-speed  constraints  with  the  same 
level  of  automation  that  previously  exlr.tsd  for  metallic  construction. 

j? ' A deflection  constrsint  generalized  as  a linear  combination  of 
translational  displacements  at  several  structural  nodea  can  be 
applied  to  metallic  and  compocite  structures.  This  permits  the 
treatment  of  constraints  on  lifting-surface  twlat  nnd  camber. 

V 3 Multiple  deflection  constraints  can  be  treated  in  a limited,  but 
practically  important,  class  of  problems  by  multiple  submissions 
of  single-constraint  cases.  This  process  has  been  largely  automated. 

"»  f Elastic  restraint  st  structural  support  points  can  be  treated. 

>o.  r Reactions  at  rigid  and  elastic  supports  are  computed  and  their 
resultants  determined,  

This  report  is  presented  in  two'yolumes : 
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Volume  I - Theory  and  Application. 

Volume  II  - Program  User's  Manual. 
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FOREWORD 

This  final  report  was  prepared  by  the  Structural  Mechanics  Section  of 
the  Grumman  Aerospace  Corporation,  Bethpage,  New  York,  for  the  Structural 
Mechanics  Division,  Air  Force  Flight  Dynamics  Laboratory,  Wright- Patterson 
Air  Force  Base,  Ohio.  The  work  was  performed  under  Contract  No.  F33615- 
77-C-3027,  which  was  initiated  under  Project  No.  2401,  "Structural  Mechanics", 
Task  No.  02,  "Design  and  Analysis  Methods  for  Aerospace  Vehicle  Structures". 
Initially  Mr.  R.  F.  Taylor  (FBR)  was  the  AFFDL  engineer  for  this  contract, 
after  which  Mr.  T.  M.  Harris  (FBR)  assumed  the  position. 

The  report  consists  of  two  volumes.  Volume  1,  entitled  "Theory  and 
Application",  describes  the  analysis  and  redesign  procedures  used  in  a com- 
puter program  system  (FASTOP-3)  for  minimum-weight  design  of  metallic  and 
composite  lifting- surface  structures  subjected  to  combined  strength  and  flutter- 
spesd  requirements  or  combined  strength  and  deflection  requirements.  Em- 
phasis is  placed  on  capabilities  that  are  new  in  FASTOP-3  compared  with 
earlier  versions  of  the  FASTOP  system.  Detailed  instructions  on  the  use  of 
the  FASTOP-3  system  are  provided  in  Volume  II,  entitled  "Program  User's 
Manual".  The  report,  which  covers  work  conducted  between  1 May  1977  and 
1 May  1978,  was  submitted  to  the  Air  Force  in  May  1978. 

Drs.  J.  Markowitz  and  G.  Isakson  were  the  Project  Engineers.  The 
authors  gratefully  acknowledge  the  suggestions  and  advice  of  Messrs.  E. 
Lerner  and  J.  Smedfjeld  and  Dr.  W.  Lansing.  They  are  also  indebted  to 
Messrs.  D.  George  and  P.  Stylianos  for  assistance  on  programming  aspects, 
and  to  Mr.  G.  Schriro  for  assistance  with  report  preparation. 
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Section  1 
INTRODUCTION 

In  December  1975,  Grumman,  under  contract  to  the  Air  Force  Flight  Dy- 
namics Laboratory,  completed  development  of  a large  scale  computer  program 
system  named  FASTOP  (Flutter  And  STrength  Optimization  Program)  which  is 
capable  of  performing  integrated  interdisciplinary  analysis  and  redesign  for 
aircraft  lifting-surface  structures  (Reference  1).  That  program,  henceforth 
referred  to  as  FASTOP- 1,  can  analyze  designs  with  respect  to  loads,  strength, 
vibration  and  flutter,  but  most  importantly,  the  program  is  also  capable  of 
obtaining  near-optimum  designs  for  lifting-surface  structures  subjected  to 
both  strength  and  flutter-speed  requirements. 

The  strength  analysis  and  optimization  module  of  FASTOP- 1 is  based, 
with  some  modifications,  on  a version  of  ASCP  (Automated  Structural  Optimiza- 
tion Program)  that  had  been  developed  for  the  Air  Force  several  years  earlier 
(Reference  2).  That  version,  now  known  as  ASOP-1,  is  limited  to  noncom- 
posite materials.  Accordingly,  FASTOP- 1 cannot  ba  used  for  interactive 
strength/flutter  resizing  of  composite  structures;  it  can  be  used  for  flutter 
resizing  alone,  but  the  associated  input  data  requirements  are  rather  cumber- 
some. 

The  ASOP-1  program  was  followed  by  the  development  of  two  successively 
more  sophisticated  versions,  ASOP-2  and  ASOP-3,  each  of  which  have  increased 
capability  for  treating  composite  materials  (References  3 and  4).  Among  other 
improvements,  ASOP-2  is  able  to  perform  strength  resizing  of  laminated  fila- 
mentary composites  of  0°/90°/i45°  balanced  layup.  The  most  recent  version, 
ASOP-3,  treats  composites  in  a rather  comprehensive  manner.  Specifically,  it 
can  accommodate  laminates  with  up  to  six  arbitrary  fiber  directions,  enforce 
"balanced"  laminate  requirements  if  desired,  and  apply  practical  design  criteria 
during  resizing.  Moreover,  the  program  has  an  improved  algorithm  for  inter- 
active strength/deflection  resizing  which  follows  closely  the  procedure  for 
strength/  flutter  resizing  in  FASTOP-1. 

Ir.  view  of  the  ASOP-3  program's  extensive  capabilities  with  respect  to 
the  treatment  of  composites,  an  effort  was  undertaken  to  introduce  those  cap- 
abilities into  the  FASTOP  system.  Two  major  tasks  were  accomplished. 
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First,  all  of  the  features  of  ASOP-3  were  incorporated,  with  no  loss  of  capa- 
bility, into  the  strength  analysis  and  optimization  module  of  FASTOP-1. 
Second,  the  flutter  optimization  module  cf  FASTOP-1,  as  v ’1  as  all  of  the 
subroutines  associated  with  interactive  strength/fluttrr  res-  , were  modified 
to  accommodate  compos  «.•  i.  The  new  version  of  FASTOP,  designated 
FASTOP-3,  is  described  n the  present  report.* 

Several  useful  features,  not  present  in  either  FASTOP-1  or  ASOP-3, 
have  also  been  introduced  into  FASTOP-3.  These  include  the  computation  of 
reactions  at  rigid  and  flexible  support  points,  and  an  automated  restart  capa- 
bility for  deflection  resizing  which  greatly  facilitates  the  treatment  of  multiple 
deflection  constraints  by  means  of  successive  applications  of  the  program. 

This  report  has  been  written  as  a supplement  to  the  original  FASTOP 
report  (Reference  1)  and  is  not  intended  to  serve  as  z stand-alone  document. 
Where  FASTOP  has  changed  little,  or  not  at  all,  the  program  description  is 
kept  very  brief,  and  it  is  expected  that  a reader  new  to  FASTOP  will  read 
the  pertinent  portions  of  Reference  1.  Where  extensive  changes  have  been 
made,  as  in  resizing  for  strength  and  deflection  constraints,  the  description 
is  sufficiently  complete  that  the  reader  can  gain  a good  understanding  of  the 
principal  features  of  the  process.  However,  for  some  details  on  strength  and 
deflection  resizing  the  interested  reader  should  refer  to  Reference  4. 

The  major  unique  capability  of  FASTOP-3,  namely  interactive  strength/ 
flutter  resizing  of  composite  lifting-surface  structures,  has  been  demonstrated 
on  two  sample  structures.  Results  of  the  two  redesign  studies  are  presented 
in  some  detail.  Also,  in  order  to  verify  that  the  ASOP-3  capabilities  have 
been  incorporated  successfully  into  FASTOP-3,  it  was  necessary  to  apply 
FASTOP-3  to  a strength/deflection  redesign  problem  that  had  been  addressed 
previously  by  ASOP-3.  The  results  obtained  for  this  additional  redesign  pro- 
blem, using  FASTOP-3,  are  included  in  this  report. 


*The  designation  FA3TOP-2  pertains  to  a modified  version  of  FASTOP-1 
developed  by  the  University  of  Dayton  and  used  by  personnel  at  the  Air 
Force  Flight  Dynamics  Laboratory. 
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Section  2 

MODIFICATIONS  TO  THE  ORIGINAL  FASTOP  SYSTEM 

This  section  contains  descriptions  of  the  theory  and  procedures  assoc- 
iated with  the  analysis  and  redesign  capabilities  of  FASTOP>3.  The  focus  is 
on  the  differences  between  FASTOP-3  and  the  original  version  of  the  program. 
Where  these  differences  are  small,  or  nonexistent,  the  discussion  is  very 
brief  and  the  reader  is  referred  to  the  proper  section  of  the  original  FASTOP 
theory  report  (Reference  1).  Where  the  differences  are  extensive,  namely  in 
the  area  of  strength/deflection  analysis  and  redesign , the  discussion  is  quite 
thorough  and  it  is  hoped  that  the  material  is  sufficient  for  the  reader  to 
grasp  the  important  concepts  and  procedures.  For  related  details,  however, 
the  reader  should  consult  the  ASOP-3  final  report  (Reference  4). 

The  following  subsections  are  arranged  so  as  to  correspond  with  the 
order  in  which  topics  are  presented  in  Reference  1 . 

2.1  OVERVIEW  OF  THE  FASTOP  SYSTEM 

The  FASTOP  system  is  divided  into  two  major  programs,  the  first  of 
which  addresses  static  analysis  and  redesign  functions,  and  the  second  of 
which  addresses  aspects  of  dynamic  analysis  and  redesign.  Each  of  these  two 
programs,  namely  the  Strength  Optimization  Program  (SOP)  and  the  Flutter 
Optimization  Program  (FOP),  is  further  divided  into  a number  of  special 
purpose  modules.  For  FASTOP-3,  the  capabilities  of  some  of  these  modules 
have  changed.  Accordingly,  the  new  capabilities  of  all  the  various  analysis 
and  redesign  modules  are  summarized  below: 

Applied  Loads 

• Aerodynamic 

- Maximum  number  of  flight  conditions  20 

- Maximum  number  of  aerodynamic  panels 

o Subsonic  150 

o Supersonic  100 

- Maximum  number  of  control  surfaces  8 
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• Inertial 


Maximum  number  of  flight  conditions 
Maximum  number  of  distributed  (point)  masses 
Maximum  number  of  concentrated  (mass  and  inertia) 
masses 


20 

1G00 

1Q0 


Strength  and  Deflection  Analysis  and  Resizing 


Analysis 

- Maximum  allowable  number  of  finite  elements  to 
define  the  structures  model 

- Maximum  number  of  struetures-model  node  points 

- Maximum  number  of  structures-model  degrees  of 
freedom 

- Maximum  number  of  applied  load  conditions 
Resizing 


- Maximum  number  of  variables  for  resizing 
o strength 
o deflection 


Vibration  Analysis 

• Applicable  to  cantilever  or  free-free  structures 


- Maximum  number  of  dynamics-model  degrees  of 
freedom 


Flutter  Analysis 

- Assumed-pressure-function  and  doublet- 
lattice  routines  for  subsonic  flow 

- Mach-box  routine  for  supersonic  flow 

- Maximum  number  of  modes  for  flutter  analysis 

- Maximum  number  of  control  surfaces  on  main 
surface  for  doublet-lattice  and  Mach-box  routines 

- Maximum  number  of  aerodynamic  panels 
o doublet-  lattice 

o Mach-box 


3000 


10C0 

6000* 


20 


10000 

6000 


200 


M=0  -0.9 


M=1 .2  -3.0 
20 
5 


400 

350 


♦Reduces  to  3000  if  su  ^sequent  flutter  resizing  uses  a free-free  vibration 
model;  unchanged  for  cantilever  model. 
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Flutter  Resizing 

- Maximum  number  of  variables  which  can  be  resized  for 
flutter 

o structural  variables  6000* 

o mass-balance  variables  20 

A comprehensive  review  of  the  FASTOP  system  is  contained  in  Section  2 of 
Reference  1 . 

2.2  APPLIED  LOADS  ANALYSIS 

The  computation  of  applied  loads  in  FASTOP-3  differs  from  that  of  the 
original  version  only  in  that  the  maximum  number  of  load  conditions  has  been 
increased  to  20.  Both  aerodynamic  and  inertial  loads  can  be  generated,  each 
type  being  computed  with  respect  to  its  own  math  model.  The  two  sets  of 
loads  are  then  transformed  to  the  structures  model  (by  means  of  transforma- 
tion matrices  discussed  elsewhere)  for  use  in  structural  analysis.  Computation 
of  the  aerodynamic-influence-coefficient  matrices  is  based  on  either  subsonic 
vortex-lattice  theory  or  supersonic  source-distribution  theory.  For  the 
latter,  two  general  flow  conditions  are  considered,  namely  supersonic  leading 
and  trailing  edges,  and  subsonic  leading  and  supersonic  trailing  edges. 
Details  of  the  applied  loads  analysis  are  presented  in  Section  3 of  Reference  1. 

2.3  STRENGTH  AND  DEFLECTION  ANALYSIS  AND  RESIZING 
2.3.1  Summary 

The  strength  and  deflection  analysis  and  resizing  capability  of  the  ASOP-3 
program  (Reference  4)  has  been  incorporated  into  the  SOP  program,  replacing 
the  earlier  version  of  ASOP  previously  used  in  FASTOP.  This  change  greatly 
enhances  the  utility  of  the  SOP  program,  as  follows: 


♦Reduces  to  4800  if  maximum  allowable  sizes  have  been  prescribed  'or  any  of 
the  structural  variables  involved  in  the  flutter  resizing  process. 
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(1)  Composite  materials  of  laminated  filamentary  construction  cen  be 
treated  explicitly  and  efficiently.  In  any  composite  element,  the 
number  of  laminae  in  each  group  of  laminae  of  a given  material,  and 
having  a common  fiber  orientation,  is  treated  as  a separate  design 
variable . 

(2)  The  program  can  resize  to  satisfy  strength  constraints  and  a single 
deflection  constraint  interactive’y . The  deflection  that  is  to  be 
constrained  can  be  generalized  in  the  sense  that  it  can  be  expressed 
as  a linear  combination  of  translational  displacements  at  selected  . 
nodes.  However,  the  constraint  con 'be  applied  in  only  a single 
load  condition.  Multiple  deflection  constraints  can  be  treated  in  a 
limited  class  of  problems  by  multiple  passes  through  the  program, 
within  each  of  which  a single  deflection  constraint  is  applied. 

(3)  Up  to  twenty  load  conditions  can  be  accommodated  in  resizing  for 
strength . 

2.3.2  Finite-Element  Modeling 

The  available  finite  elements  are  the  same  as  in  the  ASOP-3  program. 
They  consist  of  a uniform-cross-section  bar  element  (type  1),  two  prismatic 
beam  elements  (types  2 and  11),  a triangular  membrane  element  (type  4),  two 
quadrilateral  membrane  elements  (types  5 and  8),  end  a quadrilateral  shear 
panel  (type  6).  In  contrast  with  the  original  version  of  FASTOP,  all  of  the 
membrane  elements  (types  4,  5,  and  8)  can  be  of  composite  construction,  and 
the  type-2  uniform  beam  element  has  been  generalized  to  provide  for  offset  of 
the  element  nodes  from  the  beam  centroidal  axis.  The  plate  bending  elements 
(types  15  and  16)  and  combined  membrane  and  plate  bending  elements  (types 
17  and  18)  are  not  available  in  the  new  program. 

Composite  elements  can  have  up  to  six  layers,  where  a ‘'layer”  is  defined 
as  the  aggregate  of  all  laminae  of  a given  composite  material  (for  example, 
graphite/epoxy),  and  with  fibers  in  a given  direction.  The  fiber  directions 
for  the  different  layers  can  be  arbitrary,  except  that,  in  segregate,  the 
layers  should  constitute  a laminate  whose  strength  is  filament-controlled  (as 
defined  in  Subsection  2. 3. 4. 3).  It  should  be  noted  that  the  laminae  in  a 
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layer,  as  defined  above,  will  not  norraal'y  be  contiguous  in  an  actual  laminate , 
but,  in  membrane  action,  their  combined  effect  will  be  the  same  as  if  they 
were  contiguous . 

A composite  material,  as  defined  in  the  program  for  the  purpose  of  data 
input,  consists  of  the  whole  laminate,  including  numbers  of  layers,  filament 
direction  in  each  Isyer  (relative  to  the  xp-axia  in  the  property  axis  system), 
and  properties  of  the  material  in  each  layer.  Starting  and,  if  desired,  mini- 
mum and  maximum  numbers  of  laminae  in  each  layer  are  also  specified. 

Layers  of  composite  elements  are  treated  internally  in  the  program  as 
separate  elements,  except  in  resizing  for  stress  constraints,  where  the  inter- 
action between  them  is  accounted  for,  as  discussed  in  Subsection  11. 3, 4. 3. 
The  layers  of  a composite  element  thus  constitute  a stack  of  elements  that  are 
connected  at  the  comer  nodes  only.  In  the  case  of  quadrilateral  elements, 
which  consist  of  an  assemblage  of  triangular  elements  with  a common  interior 
node,  there  may  be  some  relative  displacement  of  these  interior  nodes  among 
layers  in  a laminate.  This  effect  is  believed  to  be  small,  however,  and  should 
not  introduce  any  significant  error. 

In  the  specification  of  material  properties,  three  different  classes  of 
materials  are  distinguished:  isotropic,  orthotropic,  and  composite,  with 

composite  materials  defined  as  discussed  above.  If  the  elastic  properties  of  a 
member  are  more  generally  anistropic,  they  may  still  be  introduced,  but  the 
necessary  data  must  then  be  entered  through  member  data  cards.  (The 
members  are  the  finite  elements  of  the  structural  model,  and  the  terms  "member" 
and  "finite  element"  or  ''element"  are  used  here  interchangeably.)  The  member 
data  specify  the  element  type,  the  material  code,  the  nodes  it  connects,  and 
other  required  geometrical  data  and  program  clues.  In  addition,  material 
stiffness  and  strength  data  may  be  specified,  if  desired,  with  the  member 
data,  rather  than  through  the  material  property  data,  as  explained  in  Volume 
II  of  this  report.  Likewise,  in  the  case  of  composite  members,  the  layer 
properties  may  be  introduced  on  member  data  cards,  instead  of  through  the 
material  properties  input. 

For  one-dimensional  members  or  elements  (bars  and  beams),  a single 
system  of  orthogonal  coordinate  axes  suffices.  In  the  case  of  two-dimensional 
elements,  however,  at  least  two  systems  of  orthogonal  axes  are  used,  as 
shown  in  Figure  2.1  for  a quadrilateral  element.  One  of  these,  named  the 
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"local  element  exes,"  , y£,  *£,  le  oriented  «o  that  the  origin  ia  at  node  i, 
where  the  nodes  are  entered  in  the  program  in  the  order  i,  j,  k,  1,  and  are 
arranged  as  shown  in  Figure  2.1.  The  x£-axis  is  along  the  edge  i*j,  positive 
toward  j;  the  y£-axis  is  positive  toward  the  side  on  which  the  element  lies; 
and  the  z£*axis  completes  a right-handed  triad.  While  the  direction  i to  j is 
shown  as  being  counterclockwise  around  the  element  in  Figure  2.1,  it  can  Just 
as  well  be  clockwise,  as  long  as  node  k is  on  a common  edge  with  node  i,  and 
node  1 with  node  j.  For  planar  elements,  the  local  element  axes  are  in  the 
plane  of  the  element,  while  in  the  case  of  warped  quadrilateral  elements,  they 
are  in  a plane  defined  by  a pair  of  straight  lines  joining  the  midpoints  of 
opposite  sides,  as  discussed  in  Appendix  E.2  of  Reference  4. 

The  use  of  a second  set  of  axes,  the  "property  axes,"  xp,  yp.  is  manda- 
tory when  the  material  is  not  isotropic.  These  axes,  which  are  in  the  same 
plane  as  the  x£  and  y£  axes,  are  aligned  with  directions  having  significance 
in  the  definition  of  material  properties.  For  example,  in  the  case  of  ortho- 
tropic materials,  they  are  the  axes  of  symmetry  in  the  material  properties. 


vv  vr 


//  / 

i / 

/ 


Figure  2.1  M«mb#r  Axis  SyiUms 


For  composite  materials,  they  can  be  arbitrary,  hut  are  normally  related  in 

some  simple  manner  to  the  fiber  directions,  as,  for  example,  in  the  case  of  a 

0°/90°/i45°  laminate,  where  a natural  choice  for  the  x-axis  would  be  the  0° 

P 

direction.  It  is  necessary  to  define  property  axes,  even  for  isotropic  ele- 
ments, when  the  user  desires  that  the  stress  and  strain  output  rot  be  re- 
ferred to  the  local  element  axes.  The  property  axes  then  define  the  axes  to 
which  the  stresses  and  strains  are  to  be  referred. 

The  angle  between  the  x#  and  x axes  is  referred  to  as  the  p-ungle," 

X p 

and  is  positive  as  shown  in  Figure  2.1,  that  is,  when  a rotation  from  x£  to 
Xp  is  away  from  the  element,  rather  than  into  the  element.  There  are  two 
ways  of  specifying  this  angle  for  each  member.  One  way  is  to  enter  it  directly 
with  the  data  for  each  member,  but  this  requires  that  it  be  precalculated  for 
all  the  members  - a task  that  can  be  formidable.  The  other  way  is  to  separate 
the  members  into  groups  or  "zones",  and  have  the  program  automatically 
compute  the  p-angles  for  all  the  members  in  each  zone  on  the  basis  of  a 
"reference  direction"  defined  for  that  2one.  The  0-angle  is  the  angle  between 
the  projection  of  the  reference  direction  in  the  plane  of  the  member  (or  the 
i-j-k  plane  in  the  case  of  warped  elements)  and  the  x^-axis.  The  relations 
used  >n  that  calculation  are  presented  in  Appendix  C of  Reference  4.  Only 
the  zone  number  need  then  be  specified  for  each  member,  the  reference 
directiots  being  defined  elsewhere  for  each  zone. 

In  the  case  of  composite  members,  there  is  a third  set  of  orthogonal  axes 
defined  for  each  layer.  These,  axes  are  the  "fiber  axes,"  x^,  y{,  which  are 
in  the  same  plane  as  the  Xp  and  yp  axes  (and  also  the  x^  and  yg  axes)  and 
are  aligned  so  that  the  x^-axis  is  in  the  fiber  direction.  The  angle  between 
the  x^  and  Xp  axes  is  the  angle  + defined  for  each  layer  in  the  input  for  the 
composite  material.  It  is  positive  in  the  Same  sense  as  p,  that  is,  when  a 
rotation  from  Xp  to  x^  is  away  from  the  element,  as  shown  in  Figure  2.1. 

While  the  various  systems  of  axes  have  been  defined  above  with  reference 
to  quadrilateral  elements,  the  definitions  remain  the  3ame  for  triangular  elements 

2.3.3  Analysis  Procedure 

The  analysis  procedure  is  essentially  the  same  as  in  the  original  version 
of  FASTOP.  That  is,  it  applies  the  matrix  displacement  method,  sta''';ing  the 
element  stiffness  matrices  to  form  a structure  stiffness  matrix,  with  due 


account  taken  of  boundary  conditions.  The  matrix  equation  relating  applied 
loads  to  nodal  deflections,  through  the  structure  stiffness  matrix,  is  then 
solved  for  the  nodal  deflections,  which  are  subsequently  used  to  determine 
corner  forces  acting  on  the  individual  elements. 

As  discussed  in  Subsection  2.3.2,  the  layers  of  a composite  element  are 
treated  as  separate  elements  in  the  stacking  process.  In  determining  their 
stiffness  matrices  for  this  purpose,  the  full  stiffness  properties  of  the  lamina 
material,  including  the  matrix  material,  are  taken  into  account.  This  differs 
from  the  procedure  in  the  resizing  process,  where  the  stiffness  of  the  matrix 
material  is  neglected  in  determining  layer  stresses. 

A new  capability  introduced  into  the  program  involves  the  calculation  of 
the  support  reactions.  The  force  and  moment  resultants  of  all  the  reactions 
for  each  load  condition  are  determined,  and  can  then  be  compared  with  the 
corresponding  resultants  of  the  applied  loads.  This  provides  a check  of  the 
degree  to  which  equilibrium  conditions  are  satisfied,  and,  consequently,  a 
check  of  the  numerical  precision  achieved  in  the  analysis. 

The  load-displacement  relations  may  be  expressed  in  the  matrix  form, 
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where  {6}  is  the  vector  of  unconstrained  nodal  deflections,  {Pi  is  the  vector 
of  applied  loads  in  corresponding  deg. ees  of  freedom,  {R}  is  the  vector  of 
support  reactions,  {0}  is  a null  vector,  and  the  stiffness  matrizt  is  partitioned 
in  a consistent  manner.  Thus, 

{PJ  * IKU]  {6}  (2.2) 

{R}  * IK21]  {«}  (2.3) 

Equation  (2.2  if.  solved  for  {6},  which  is  then  substituted  into  Equation 
(2.3)  to  yielo  the  support  reactions.  As  in  the  original  version  of  FASTOP, 
Equation  (2.2)  is  solved  by  employing  the  ChoJesky  algorithm  for  decomposition 
of  positive-definite  symmetric  matrices  (see  Reference  1). 
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2.3.4  The  Stress-Constraint  Mode 


2.3.4. 1 Basic  Procedure  and  Eiement  Stress  Determination 

Two  different  phases,  or  "modes’* , in  the  redesign  process  are  distin- 
guished. In  the  "stress-constraint  mode,"  a number  of  cycles  of  analysis  and 
resixing  for  ttvess  and  minimum-  and/or  maximum-gage  constraints  are  perform- 
ed, until  > ,«vergence  criterion  is  satisfied  or  the  number  of  cycles  has 
reached  a specified  maximum.  The  design  should  then  be  fully  stressed  or 
nearly  fully  stressed.  If  deflection-constraint  resizing  is  not  to  be  done,  the 
process  is  completed.  If  deflection-constraint  resizing  is  to  be  done,  there  is 
a subsequent  mode,  referred  to  as  the  "deflection-constraint  mode,"  in  which 
stn.ss-constraint  resizing  and  deflection-constraint  resizing  are  done  inter- 
actively, and  in  which  minimum-  and/or  maximum-gage  constraints  may  also  be 
applied.  It  should  be  noted  that  minimum  gages  usually  represent  limitations 
associated  with  practical  construction,  while  the  maximum  gage  is  normally  used 
as  a means  of  fixing  an  element's  gage  when  the  minimum  and  maximum  gages 
are  set  equal  to  one  another. 

The  basic  procedure  for  resizing  based  on  allowable  stresses  and  member 
gage  limitations  is  essentially  the  same  as  in  the  original  version  of  FASTOP. 
The  stiffness  matrix  of  each  finite  element  in  the  structural  model  is  assumed 
to  be  linearly  related  to  a single  design  variable  for  that  element  (bar  cross- 
sectional  area,  skin  gage,  etc.).  In  the  case  of  the  beam  elements,  it  is 
assumed  that  the  radii  of  gyration  of  the  cross  section  remain  constant  during 
resizing,  so  that  the  moments  of  inertia  of  the  cross  section,  and  conse- 
quently the  bending  stiffness,  are  proportional  to  the  cross-secticnal  area, 
which  is  the  design  variable  in  this  case. 

A:  initial  design,  specifying  values  of  the  design  variables,  is  selected, 
and  an  analysis  to  determine  nodal  deflections  and  element  stresses  is  carried 
out,  for  a given  set  of  applied  loading  conditions.  The  state  of  stress  in 
each  element,  for  each  loading  condition,  is  then  used  in  conjunction  with  a 
failure  criterion  to  determine  a "stress  ratio" . This  ratio  provides  a measure 
of  the  extent  to  which  the  stress  constraint  is  satisfied  or  violated,  and  is 
discussed  in  detail  in  Subsections  2. 3. 4. 2 and  2. 3. 4. 3.  It  is  equal  to  1.0  if 
the  failure  criterion  is  satisfied  exactly.  The  maximum  value  of  the  stress 
ratio,  for  all  loading  conditions  for  each  element,  is  then  used  as  a multi- 


plying  factor  In  resizing  the  deiign  variable  for  that  element.  The  procedure 
is  considerably  more  complicated  in  the  case  of  composite  elements,  as  ex* 
plained  in  Subsection  2. 3. 4. 3. 

This  process  of  analysis  and  redesign  is  repeated  cyclically,  until  a 
specified  number  of  cycles  have  been  performed,  or  until  a convergence  cri- 
terion is  satisfied.  A converged  design  is  referred  to  as  a fully  stressed 
design;  that  is,  a design  in  which  each  element  is  either  stressed  to  the 
maximum  allowable  extent  in  at  least  one  loading  condition  (without  being 
overstressed  in  any  loading  condition),  or  is  at  a minimum  or  maximum  pre- 
scribed gage. 

The  nodal  stress  method  for  element  stress  determination,  which  was 
used  in  the  original  version  of  FASTOP,  has  been  found  to  have  shortcomings, 
and  has  been  replaced  in  FASTOP-3.  Average  stresses  are  now  determined 
for  an  element  directly  from  the  nodal  forces  acting  on  that  element.  In  the 
case  of  the  bar  element,  which  is  a uniform-strain  element,  the  stress  is 
simply  the  quotient  of  the  nodal  force  (which  is  the  uniform  axial  toad  in  the 
bar)  and  the  cross-sectional  area  of  the  bar.  The  triangular  membrane 
element  is  similarly  based  on  the  assumption  of  uniform  strain,  and  the  aver- 
age stress  in  it  is  simply  the  uniform  stress  associated  with  the  uniform 
strain.  The  matrix  transformation  relation  between  average  stress  and  the 
comer  forces  in  that  element  is  derived  in  Appendix  E of  Reference  4. 

The  determination  of  average  stress  in  the  quadrilateral  membrane  ele- 
ments  (types  5 and  8)  is  not  as  straightforward,  as  these  elements  are  con- 
structed of  four  uniform-strain  triangular  elements,  within  each  of  which  the 
strain  is  generally  different.  There  are  a variety  of  ways,  necessarily  ap- 
proximate in  nature,  in  which  average  stress  can  be  defined  in  such  an 
element.  The  definition  used  in  FASTOP-3  is  based  on  equilibrium  considera- 
tions. A derivation  of  the  applicable  relations  is  presented  in  Appendix  E of 
Reference  4.  The  average  shear  stress  in  the  shear  panel  (element  type  6) 
is  determined  in  the  same  way. 

In  the  beam  elements  (types  2 and  11),  the  axial  load  and  corresponding 
stress  are  determined  in  the  same  way  as  in  the  bar  element.  Rending  mo- 
ments, however,  are  determined  separately  at  each  end  of  the  element,  and. 
because  it  is  assumed  that  the  element  is  loaded  only  at  its  ends,  the  maximum 
bending  moment  occurs  at  one  of  the  ends.  The  shear  force  and  torsional 
moment,  both  uniform  along  the  element,  are  also  determined. 
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2 3. 4. 2 Resizing  Algorithm  for  Noncomposite  Elements 

the  resizing  of  bar  elements  la  particularly  simple  because  of  the  uniaxial 
stress  state  in  tlu-m,  the  stress  ratio  being  the  ratio  of  the  actual  stress  to 
the  allowable  stress.  In  the  resizing  of  beam  elements,  it  is  assumed  that  the 
element  is  loaded  in  bending  primarily  about  only  one  of  the  two  transverse 
axes:  the  z-axis.  The  bending  moment  about  that  axis  is  then  determined  at 
the  two  ends  of  the  -dement,  and  corresponding  extreme-fiber  stresses  are 
determined  at  the  two  ends,  assuming  that  the  distance  of  the  extreme  fiber 
from  the  neutral  axis  is  equal  to  the  radius  of  gyrathn  of  the  cross  section 
(which,  in  effect,  asst  tees  that  the  bending  material  ie  concentrated  at  the 
extreme  fiber).  The  beading  stresses  are  then  combined  with  the  stress  due 
to  the  axial  load,  yielding  four  values  of  stress  - two  extreme-fiber  stresses 
at  each  end  of  the  element  - and  corresponding  stress  ratios  are  determined. 
The  largest  of  these  stress  ratios  is  then  selected  for  resizing  purposes. 

The  biaxial  stress  state  in  membrane  elements  requires  that  a failure  cri- 
terion, providing  for  the  interaction  between  stress  components,  be  used.  In 
the  case  of  isotropic  materials,  it  is  common  to  use  the  von  Mises  yield  cri- 
terion, with  ultimate  allowable  stresses  usually  replacing  the  yield  stresses, 
and  that  criterion  is  used  in  FASTOP-3,  in  the  following  modified  form: 


» 1 (2.4) 


where  Fx  is  either  the  allowable  tensile  stress  Ft  or  the  allowable  compressive 
stress  Fc,  depending  upon  the  sign  of  ox;  F^  is  either  F{  or  Fc>  depending 
upon  the  sign  of  a ; and  F is  the  allowable  shear  stress.  Ft,  F and  F are 
always  taken  as  positive  quantities.  The  stress  ratio  in  this  case  is  the  left- 
hand  side  of  Equation  (2.4). 

In  the  case  of  orthotropic  materials,  the  picture  is  less  clear  and  there 
is  no  universally  accepted  failure  criterion.  Two  relatively  simple  criteria, 
having  a somewhat  rational  basis,  are  Hill's  generalization  of  the  von  Mises 
criterion  and  a criterion  developed  originally  at  the  Forest  Products  Labora- 
tory (References  5 and  6).  In  the  absence  of  conclusive  experimental  evi- 
dence favoring  either  one  of  these  criteria,  the  latter  was  selected  for  use  in 
ASCP-3,  and  is  used  also  in  FASTOP-3. 
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where  a*.  oy  and  are  the  atreaa  components  in  the  pivperty  axis  system; 

F , F and  F are  the  corresponding  allowable  stresses  in  the  absence  of  the 
* y * 

other  components;  and  failure  is  presumed  to  occur  when  any  one  of  the  three 

relations  is  satisfied.  F and  F are  tensile  or  compressive  allowable 

* y 

stresses,  as  appropriate,  and  are  always  taken  as  positive  quantities. 

In  applying  this  crit'rton,  the  left-hand  sides  of  all  three  relations 
are  evaluated,  and  the  largest  of  the  three  becomes  the  governing  stress 
ratio,  which  is  then  used  as  a multiplying  factor  in  resizing,  as  described 
in  Subsection  2.3.4  l. 

The  failure  criterion  expressed  by  Equation  (2.5)  is  shown  graphically 

as  an  envelope,  in  Figure  2.2,  for  the  case,  i » 0.  The  first  of  the  equa- 

xy 

tions  is  represented  by  the  ellipse,  and  the  remaining  two  equations  effectively 
apply  cutoffs  in  the  first  and  third  quadrants.  It  can  be  seen  that  the 
ellipse,  in  Figure  2.2,  also  represents  Equation  (2.4),  which  ts  the  failure 
criterion  for  isotropic  materials.  If  desired,  the  cutoffs  shown  in  Figure  2.2 
can  be  applied  to  isotropic  materials  as  well,  by  suitable  .adjustment  of  the 
input  data,  a*  explained  in  Volume  II  of  this  report. 

The  new  element  gages  obtained  in  this  way  are  checked  against  the 
minimum  gages  specified  in  the  input  data,  and  are  revised  accordingly,  if 
necessary.  If  the  stress  resizing  is  being  done  as  part  of  a cycle  in  the 
deflection-constraint  mode,  as  discussed  in  Subsection  2. 3. 4. 2,  a similar 
check  is  made  against  the  minimum  gages  established  by  deflection-constraint 
resizing,  and  necessary  revisions  are  made.  Similarly,  if  the  stress  resizing 
is  being  done  subsequent  to  flutter  resizing,  the  flutter-governed  gages  are 
now  treated  as  minimum  gages.  A check  is  then  made  against  the  maximum 
gages  specified  in  the  input  data,  and  necessary  revisions  are  made. 
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Flfur*  2.2  F«lkw«  Criterion  for  Ortho  tropic  Meter  Ml  in  Biaxial  Straw 

While  composite  elements  can  (it  is  expected  that  they  usually  will)  be 
resized  by  the  comprehensive  ptocedure  described  in  the  following  section,  a 
cruder  approach  can  be  used,  i’.t  which  they  are  treated  as  orthotropic  elements 
and  resized  by  the  procedure  described  above.  In  that  case,  it  is  necessary 
that  the  layup  be  such  that  the  relative  numbers  of  laminae  in  the  various 
fiber  directions  are  maintained  constant,  so  that  appropriate  allowable  stresses 
can  be  selected.  It  will  be  seen  that  the  cutoffs  described  in  the  following 
section  bear  some  similarity  to  those  shown  in  Figure  2.2. 

2. 3. 4. 3 Resizing  Algorithm  for  Composite  Elements 

The  criteria  governing  the  failure  of  composites  are  more  complex  than 
those  governing  the  failure  of  noncomposite  materials.  In  consequence,  the 
algorithm  for  composite  element  resizing  is  necessarily  more  complex,  requiring 
that  the  laminate  be  treated  as  a unit,  so  that  interaction  between  layers 
may  be  properly  taken  into  account.  Furthermore,  because  of  limited  opera- 
tional experience  with  composite  materials,  it  is  desirable  to  make  certain 
conservative  assumptions  concerning  their  strength  behavior. 
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For  example,  local  cracking  or  crazing  in  the  matrix  material  may  greatly 
reduce  its  effectiveness  as  a load-carrying  agent,  even  though  it  continues  to 
serve  its  central  purpose  as  a binding  agent.  For  that  reason,  the  assumption 
is  made  here  that  all  the  load  is  carried  by  the  filaments.*  This  assumption 
ts  applicable  only  to  so-called  filament-controlled  composites:  those  in  which 

the  layup  is  such  that  filament  directions  are  sufficient  in  number  and  distri- 
bution that  any  component  of  the  laminate  stress  resultant  can  be  resisted  by 
filaments  alone  It  would  not  be  a tenable  assumption  in  the  case  of  matrix- 
controlled  composites:  those  in  which  the  load-carrying  capability  of  the 

matrix  is  relied  upon.  The  program  does  not,  at  present,  accommodate  com- 
posites of  this  latter  type,  although  it  can  be  further  developed  to  do  so. 

Another  conservatism,  one  that  may  be  applied  as  an  option  in  the  pro- 
gram. relates  to  interaction  between  laminate  layers  that  is  somewhat  akin  to 
the  effect  of  hydrostatic  stress  in  metals  or  other  nominally  homogeneous  ma- 
terials It  is  known,  for  example,  that,  if  components  N and  N of  the 

x y 

laminate  stress  resultant  are  present  and  are  of  the  same  sign,  some  layers 

will  be  less  severely  stressed  than  if  either  N or  N were  absent.  Some 

x y 

designers  feel  that  it  is  unconservative  to  take  advantage  of  this  fact  and  to 
base  the  design  on  the  simultaneous  presence  of  both  components,  particularly 
since  the  prediction  of  applied  loads  is  hardly  an  exact  science,  and  the 
dynamic  nature  of  some  loading  conditions  suggests  the  possibility  that  differ- 
ent components  of  the  internal  loading  may  not  be  applied  simultaneously. 
The  program  option  referred  to  f.s  the  "cutoff  option"  makes  it  possible  to 
provide  additional  stress  checks,  with  and  successively  set  to  zero, 
and  to  use  the  results  as  additional  information  in  the  resizing  process. 

In  addition  to  filament  failure  in  tension  or  block  compression,  the  possi- 
bility of  failure  in  the  so-called  "microbuckling"  mode  may  be  taken  into 
account.  In  that  mode,  there  is  a highly-localized  buckling  of  the  filaments 
because  of  then  own  low  bending  stiffness  and  the  limited  shear  stiffness  of 
the  matrix  material,  which  is  relied  upon  to  resist  such  buckling  (Reference 


♦Note  that,  as  stated  in  Section  2.3.3,  this  assumption  is  made  only  in  the  re- 
sizing process.  In  the  analysis  of  the  whole  structure,  to  determine  nodal 
displacements  and  internal  loads,  the  full  stiffness  properties  of  the  com- 
posite elements,  including  the  contribution  of  the  matrix  material,  are  used. 
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7).  Theoretically , it  has  been  found  that  the  allowable  stress,  Gz>  in  this 

mode,  should  be  equal  to  G _/(l  - V-),  where  G_  is  the  shear  stiffness  of  the 

mi  m 

matrix  material,  and  is  the  volume  fraction  of  fibers  (Reference  7);  however, 
experimental  evidence  indicates  that  it  is  better  to  use  a value  based  on 
experimental  data  (Reference  8).  If  experimental  data  are  not  available,  the 
theoretical  value  may  be  multiplied  by  an  empirical  coefficient  which,  at  least 
for  the  case  of  boron/epoxy,  can  be  given  the  value  0.63  (Reference  9). 

To  properly  account  for  the  interaction  between  layers,  the  resizing  al- 
gorithm for  composite  elements  treats  the  entire  laminate  as  a unit  and  applies 
a convergent  iteration  procedure.  This  procedure  is  summarized  in  flow-chart 
form  in  Figure  2.3  and  is  described  in  detail  in  Reference  4.  Essentially,  it 
uses  the  stress  resultant  on  the  whole  laminate  to  determine  stresses  in  the 
layers  (neglecting  the  effect  of  matrix  stiffness),  and  then  resizes  the  layers 
using  stress  ratios.  Minimum/maximum  constraints  are  then  applied  and  the 
layer  thickness  is  rounded  up  to  the  nearest  multiple  of  the  lamina  thickness. 
The  process  is  repeated  cyclically  until  convergence  is  achieved.  A micro- 
buckling  failure  criterion,  cutoffs  (as  discussed  earlier),  and  layer  'balancing1 
may  be  applied.  Balancing  forces  specified  layers  always  to  have  the  same 
number  of  laminae.  It  is  accomplished  by  setting  the  number  of  laminae  in  all 
such  layers  equal  to  the  number  of  laminae  in  the  most  severely  stressed 
layer. 

2. 3. 4. 4 Resizing  of  Compression  Members 

The  limited  capability  introduced  into  the  original  version  of  PASTOP  for 
the  resizing  of  bars  and  shear  or  compression  panels  subject  to  buckling 
failure  has  been  retained  in  FASTOP-3.  It  involves  the  introduction  of  "sta- 
bility tables"  relating  allowable  stresses  to  internal  loading.  This  capability 
is  limited  to  noncomposite  elements . 

2.3.5  The  Deflection  Constraint  Mode 


2.3.5. 1 Resizing  Algorithm 

It  is  shown  in  Appendix  G.l  of  Reference  4 that,  for  the  case  of  a 
single  deflection  constraint  and  in  the  absence  of  other  constraints,  a minimum 


weight  is  achieved  when  the  partial  derivative  of  the  constrained  deflection 
with  respect  to  element  weight  has  the  same  value  for  all  elements,  that  is, 
when 

36 

» K (i  = 1,  2 n)  (2.6) 

3w. 

> 

where  6 is  the  deflection  to  be  constrained,  w.  is  the  weight  of  the  i1^  element 
(of  a total  of  n elements),  and  K is  a constant. 

When,  as  in  most  practical  designs,  there  are  strength  constraints  and 
minimum  or  maximum  gage  constraints  in  addition  to  the  deflection  constraint, 
the  uniform-derivative  criterion,  expressed  in  Equation  (2.6),  can  be  applied 
to  the  set  of  all  elements  not  governed  by  these  other  constraints,  the  corres- 
ponding element  weights  being  referred  to  as  the  "active"  variables.  In  that 
case,  the  criterion  is  less  rigorously  applicable,  but  should  still  give  a design 
of  nearly  minimum  weight,  as  discussed  in  Appendix  G.l  of  Reference  4. 

The  minimum-weight  design  cannot  be  arrived  at  directly.  It  is  necessary 
to  employ  an  iterative  process,  which  has  been  found  to  converge  rapidly  in 
practical  cases.  This  process  closely  resembles  the  procedure  used  in  the 
program  for  optimization  to  satisfy  a constraint  on  flutter  velocity . 

Starting  with  a nonoptimum  design,  which  may  or  may  not  satisfy  the 
prescribed  deflection  constraint,  the  following  recursion  relation,  derived  in 
Appendix  G.2  of  Reference  4,  is  applied  in  successive  cycles: 


where  w.  is  the  weight  of  the  i th  element  prior  to  resizing  in  the  current 
‘old 

cycle;  w.  is  the  weight  of  the  i th  element  following  resizing  in  the  current 
1 new 

cycle;  (36/aw.)^  is  the  partial  derivative  of  the  constrained  deflection  with 
respect  to  w.,  computed  for  the  design  existing  prior  to  resizing  in  the  current 
cycle;  and  (a6/aw)target  *s  a quantity  given  the  name  "target  derivative"  and 
is  defined  below.  The  basis  for  Equation  (2.7)  is  discussed  in  Appendix  G of 
Reference  4. 


At  the  optimum  design,  the  target  derivative  will  be  the  constant  K in 
Equation  (2.6),  and  the  derivatives  86/8w.  will  all  be  equal  to  it.  However, 
prior  to  convergence  to  an  optimum  design,  the  derivatives  86/8Wj  will  differ 
from  each  other  in  value,  and,  in  fact,  may  differ  in  sign.  Depending  upon 
the  sign  of  the  target  derivative,  some  of  these  derivatives  may  then  yield  a 
negative  value  for  the  quantity  under  the  radical  in  Equation  (2.7).  The 
corresponding  elements  will  then  have  to  be  excluded  when  Equation  (2.7)  is 
applied . 

As  the  value  of  K in  Equation  (2.6)  is  not  known  until  the  optimum  design 
is  achieved,  it  is  necessary  to  find  a value  for  the  target  derivative  that, 
when  introduced  into  Equation  (2.7),  will  yield  a design  that  satisfies  the 
constraint,  at  least  approximately.  This  is  done  by  a trial  procedure,  in 
which  a value  of  the  target  derivative  is  sought  that  will  yield  a design 
satisfying  the  relation: 

n 

^desired  = fiold 

i=l 

where  6^  is  the  value  of  the  constrained  deflection  prior  to  resizing;  desired 
is  the  desired  value  of  the  constrained  deflection  in  a resizing  step;  and  the 
summation  is  over  all  elements  of  the  model,  including  those  to  which  Equation 
(2.7)  is  not  applied. 

Equation  (2.8)  is  seen  to  provide  a second-order  approximation  (in  the 
Taylor-series  sense)  to  the  desired  value  of  the  deflection.  An  exact  value 
could  have  been  obtained  hy  a trial  procedure  in  which  a structural  analysis 
is  performed,  for  the  design  corresponding  to  each  trial  value  of  the  target 
derivative,  to  determine  the  deflection  subject  to  constraint.  However,  as 
this  operation  would  have  to  be  performed  a number  of  times,  for  successive 
trial  values  of  the  target  derivative,  and  is  expensive  computationally,  it  is 
highly  advantageous  and,  in  practice,  satisfactory  to  use  Equation  (2.8) 
instead. 


I j^8Wj)old+  (awjtargetj  ("'new'  "‘old) 
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It  was  stated  above  that  Equation  (2,7)  is  applied  to  that  group  of  ele-  1 

ments  with  derivatives,  36/3w.,  that  are  of  the  same  sign  as  the  target  deri-  > 

vative.  The  determination  of  that  sign  is  now  considered.  The  sign  is  , 

established  upon  entry  into  the  deflection-constraint  mode  and  depends  upon 
whether  the  constraint  value  of  the  subject  deflection  exceeds,  or  is  less 
than,  its  current  value  (for  the  design  existing  upon  entry  into  the  deflection- 
constraint  mode).  If  the  constraint  value  exceeds  the  current  value  alge- 
braically, and  the  constraint  is  either  (1)  an  equality  constraint  or  (2)  an 
inequality  constraint  that  has  been  violated,  then  an  increase  in  deflection  is 
desired,  and  the  proper  derivative  sign  is  that  which  is  associated  with  an 
increase  in  deflection  resulting  from  an  increase  in  element  weight,  that  is,  a 
positive  sign.  For  those  elements  with  negative  derivatives,  a reduction  in 
element  weight  will  move  the  deflection  in  the  desired  direction,  and  the 
design  variables  for  these  elements  can  be  permitted  to  decrease,  to  the 
extent  permitted  by  other  constraints.  (This  is  discussed  further  in  Sub- 
section 2. 3. 5. 2.)  In  the  reverse  situation,  where  the  constraint  value  is  less 
than  the  current  value,  Equation  (2.7)  is  applied  to  those  elements  with 
negative  derivatives.  Where  the  constraint  is  an  inequality  constraint  and  is 
not  violated  by  the  design  existing  at  exit  from  the  stress-constra jit  mode, 
r.o  further  resizing  is  necessary. 

The  sign  of  the  target  derivative  will,  by  definition,  be  the  same  as  that 
of  the  derivatives  (SS/Sw.)^  introduced  into  Equation  (2.7).  It  remains  to 
find  a value  of  the  target  derivative  that  will  satisfy  Equation  (2.8).  This  is 
done  by  taking,  as  an  initiar  trial  value,  upon  entry  into  the  deflection- 
constraint  mode,  a value  equal  to  80%  of  the  average  of  all  the  derivatives 
having  the  proper  sign,  determined  as  explained  above.  (In  subsequent  re- 
design cycles  of  the  iterative  redesign  process,  as  explained  below,  the 
starting  value  of  the  target  derivative  is  the  last  value  computed  in  the 
preceding  cycle.)  The  target  derivative  is  then  incremented  until  a value  is 
achieved  that  satisfies  Equation  (2.8),  within  a tolerance  specified  by  the 
user. 

The  value  of  6^esjrecj  in  Equation  (2.8)  is  not  necessarily  the  constraint 
value.  It  may  be  advantageous,  in  some  situations,  to  move  from  the  initial 
value  of  the  subject  deflection  to  the  vicinity  of  the  constraint  value  in  a 
series  of  shorter  steps,  rather  than  in  a single  step,  Equation  (2.8)  then 
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provides  t closer  approximation  in  each  step.  Furthermore,  when  the  con- 
straint boundary  in  the  design  apace  is  reached  by  this  procedure,  it  may  be 
at  a point  considerably  closer  to  the  optimum  design  point.  Accordingly,  the 
program  provides  an  option  that  permits  the  change  from  the  initial  value  of 
the  deflection  to  the  constraint  value  to  be  made  in  a number  of  approximately 
equal  increments,  that  number  being  selected  by  the  user. 

The  deflection  that  is  subject  to  constraint  may  be  generalised  in  the 
sense  that  it  may  be  represented  as  a linear  combination  of  nodal  displace- 
ments in  specified  degrees  of  freedom.  Thus,  for  example,  an  angular  dia- 
placement  constraint  may  be  treated  by  representing  it  as  the  difference 
between  the  translational  displacements  of  two  specified  points,  divided  by*  the 
distance  between  them.  The  two  points  specified  need  not  be  at  nodes;  their 
displacements  can  be  obtained  by  interpolation  between  nodal  displacements. 
Similarly,  a given  amount  of  camber  of  a lifting  surface,  at  a given  spanwise 
station,  can  be  specified  as  a constraint  by  a similar  representation  as  a 
linear  combination  of  nodal  displacements. 

When  composite  elements  are  included  in  the  model,  each  layer  of  such 
elements  is  treated  internally  in  the  program  as  a separate  element.  Accord- 
ingly, Equation  (2.7)  is  applied  to  individual  layers,  with  the  deflection 
derivative  being  computed  for  each  layer  with  respect  to  that  layer's  weight. 
However,  the  derivatives  of  layers  that  are  to  be  balanced  are  averaged  and 
the  average  value  is  used  in  place  of  the  actual  values.  This  approach  main- 
tains layer  balance  and  properly  accounts  for  the  effect  of  design  changes  on 
the  generalized  deflection.  In  the  trial  process  of  finding  a value  of  the 
target  derivative  that  satisfies  Equation  (2.8),  the  layer  thickness  is  not 
rounded.  Rounding  is  done  only  after  that  process  is  completed,  when  the 
layer  thickness  is  rounded  up  or  down  to  the  nearest  multiple  of  the  lamina 
thickness.  By  rounding  both  up  and  down,  the  effect  of  rounding  on  the 
constrained  deflection  can  be  minimized. 

The  evaluation  of  the  deflection  derivatives  is  carried  out  in  the  same 
way  as  in  ASOP-3,  and  is  explained  in  Appendix  G.3  of  Reference  4. 
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2.3.S.2  Resizing  Algorithm  for  interacting  Deflection,  Stress,  and  Minimum/ 
Maximum  Gage  Constraints 

As  stated  above,  the  deflection-constraint  mode  is  entered  only  if  (1) 
the  deflection  constraint  is  an  equality  constraint  (6  = desired  or  (21  it  is 
an  inequality  constraint  (^desired  or  *»®<jesired^  *s  v‘°*ate<^  b/  the 
design  existing  at  the  end  of  the  stress-constraint  mode.  (It  should  be  noted 
that  the  subject  deflection  6 and  its  desired  value  ®£jesjrecj  are  algebraic 
quantities,  not  absolute  values.)  It  is  desirable  that  the  design  existing  at 
entry  to  the  deflection-constraint  mode  should  be  very  nearly  a fully  stressed 
design . 

If  the  constraint  is  an  inequality  constraint  that  is  vio.'ated  by  the 
design  existing  at  entry  to  the  deflecjon-constraint  mode,  it  is  treated  as  an 
equality  constraint  in  the  deflection-constraint  mode.  It  is  shown  in  Appendix 
G.4  of  Reference  4 that  this  approach  should  yield  a design  that  is  near 
optimum  for  the  inequality  constraint. 

Prior  to  entry  into  the  iterative  redesign  proces.,  a determination  is 
made  of  the  algebraic  sign  of  the  derivatives  to  be  introduced  into  Equation 
(2.7),  as  discussed  in  Subsection  2.3.5. 1.  If  the  current  value  of  the  de- 
flection subject  to  constraint  is  smaller  (algebraically)  than  the  constraint 
value,  positive  derivatives  are  taken.  If  the  current  value  of  the  deflection 
subject  to  constraint  is  larger  (algebraically)  than  the  constraint  value, 
negative  derivatives  are  taken.  Once  this  determination  is  made,  it  remains 
unchanged  throughout  the  remainder  of  the  procedure,  for  reasons  that  will 
now  be  discussed. 

As  long  as  movement  from  the  initial  value  of  the  constrained  deflection 
(upon  entry  into  the  deflection-constraint  mode)  toward  the  constraint  value 
is  in  the  same  direction,  it  is  clear  that  this  sign  should  not  be  changed. 
However,  consider  the  situation  where  the  constraint  value  is  overshot,  and 
movement  in  the  reverse  direction  becomes  necessary?  If  the  sign  were  to  be 
reversed,  all  those  elements  previously  resized  by  the  deflection  constraint 
would  be  suddenly  relieved  of  such  constr.  int,  and  their  gages  could  drop  to 
values  determined  by  other  constrain' s.  Under  these  circumstances,  large 
changes  could  be  expected  to  result  from  a need  for  minor  adjustments,  as 
the  amount  of  overshoot  would  normally  be  small.  These  large  changes  could 
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be  expected  to  preclude  satisfactory  convergence  to  an  optimum  design. 
Maintaining  the  same  sign  keeps  these  adjustments  essentially  within  the  same 
group  of  elements  that  have  previously  been  governed  by  the  deflection 
constraint. 

Deflection-constraint  resizing  and  stress-constraint  resizing  are  done 
sequentially  within  each  cycle  in  the  deflection-constraint  mode,  with  an 
analysis  following  each  type  of  resizing.  There  are  thus  two  analyses  per- 
formed in  each  cycle.  This  cycling  is  continued  until  a convergence  criterion 
is  satisfied  or  the  number  of  cycles  has  reached  a specified  maximum. 

Following  stress-constraint  resizing,  the  member  gages  determined  on  the 
basis  of  the  stress  constraints  are  treated  as  minimum  gages  in  the  subsequent 
resizing  for  the  deflection  constraint.  Similarly,  following  deflection-constraint 
resizing,  the  member  gages  determined  on  the  basis  of  the  deflection  are 
treated  as  minimum  gages  in  the  subsequent  resizing  for  stress  constraints. 
Continued  to  convergence,  this  process  effectively  takes  into  account  the 
interaction  between  the  deflection  constraint  and  stress  constraints.  It  can 
be  seen  that,  upon  entry  into  the  first  cycle  of  the  deflection-constraint 
mode,  deflection-constraint  resizing  cannot  reduce  any  gages,  as  they  are 
then  all  governed  either  by  stress  constraints  or  by  minimum  gage  constraints. 
However,  in  subsequent  cycles,  reductions  can  be  made  in  gages  previously 
governed  by  the  deflection  constraint.  Such  reduction,  however,  is  not 
permitted  to  bring  a gage  below  the  largest  of  the  following:  the  gage  as 

determined  on  the  basis  of  stress  constraints,  the  specified  minimum  gage, 
and  a "MAX  CUT"  value. 

The  "MAX  CUT"  value  is  determined  by  multiplying  the  gage  prior  to  re- 
sizing by  a factor  that  is  specified  by  the  user  and  has  a . value  between  0 
and  1.  Its  intent  is  to  provide  the  user  with  a means  of  limiting  the  reduction 
that  a single  resizing  step  can  effect  in  any  gage,  in  the  event  that  excessive 
reductions  cause  an  instability  in  the  iterative  solution . 

The  algorithm  applied  in  the  deflection-constraint  mode  is  illustrated  in 
flow-chart  fora  in  Figures  2.4  and  2.5  and  is  explained  in  detail  in  Reference 
4.  Figure  2.4  shows  the  algorithm  for  alternating  deflection-constraint  and 
stress-constraint  resizing,  wh'.’e  Figure  2.5  focuses  on  the  deflection-constraint 
resizing  portion  of  the  cycle. 
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The  alternating  application  of  deflection  and  stress  constraints,  with  the 
gages  determined  by  one  being  used  as  minimum  gages  in  the  other,  as  de- 
scribed  above,  converges  to  a design  in  which  there  are  two  classes  of  ele- 
ments (or  layers  in  the  case  of  composites).  One  class  comprises  elements 
that  are  fully  stressed  or  are  at  minimum  or  maximum  specified  gage.  The 
other  class  comprises  elements  that  are  governed  by  the  deflection  constraint. 
In  this  latter  class,  the  derivatives  of  deflection  with  respect  to  element  (or 
layer)  weight  all  have  nearly  uniform  values.  Departures  from  uniformity  are 
due  to  lack  of  convergence  or,  in  the  case  of  composite  element  layers,  to 
rounding  to  an  integral  number  of  laminae.  Under  these  circumstances,  it 
can  be  expected  that  the  design  will  be  close  to  optimum,  at  least  in  a local, 
if  not  in  a global,  sense. 

It  should  be  noted  that  an  option  in  the  program  permits  the  designation 
of  selected  members  as  noncandidates  for  deflection-constraint  resizing,  al- 
though they  are  still  subject  to  the  application  of  stress  constraints  in  each 
cycle  in  the  deflection-constraint  mode,  unless  their  gages  are  being  explicitly 
fixed.  This  option  is  useful,  for  example,  when  designing  structures  having 
fixed-gage  honeycomb  core  substructure  that  is  idealized  using  shear  panels. 
Another  useful  application  is  described  in  the  following  subsection . 

2.3.S.3  Treatment  of  Multiple  Deflection  Constraints  as  a Succession  of 
Single-Constraint  Problems 

Although  only  a single  deflection  constraint,  in  a single  loading  condi- 
tion, can  be  treated  in  one  submission  of  the  program,  as  discussed  above,  it 
is  possible  to  treat  special  cases  of  multiple  deflection  constraints  by  making 
multiple  submissions  of  the  program.  The  special  cases  are  those  in  which 
the  constraints  can  be  ordered,  so  that  the  first  constraint  can  be  satisfied, 
after  which  a portion  of  the  structure  can  be  frozen  in  design  to  prevent 
further  change  in  the  corresponding  deflection;  then  the  second  constraint 
can  be  satisfied  by  redesign  of  the  remaining  structure,  after  which  a portion 
of  that  structure  can  be  frozen  in  design,  to  prevent  further  change  in  the 
corresponding  deflection;  and  so  on,  until  the  last  constraint  is  satisfied. 
Clearly,  cantilever  structures,  particularly  slender  ones,  w;:'i  deflection 
constraints  such  as  angles  of  twist  applied  at  two  or  more  stati...  along  the 
span,  fit  this  situation  to  some  degree  of  approximation. 


The  program  must  be  submitted  for  execution  as  many  times  as  there  are 
constraints  to  be  satisfied.  In  each  submission,  the  design  variables  are 
initialized  at  the  final  values  they  had  in  the  preceding  submission,  and,  in 
the  case  of  those  design  variables  that  are  nominally  to  be  frozen,  their 
newly  initialized  values  are  also  their  minimum  values.  The  word  "nominally" 
is  used  because,  while  these  design  variables  are  removed  from  candidacy  for 
deflection-constraint  resizing,  they  are  not  truly  frozen,  as  it  is  still  neces- 
sary to  apply  stress  constraints  to  them  if  overstress  is  not  to  occur.  This 
may  have  the  effect  of  further  altering  the  deflections  already  set  at  their 
constraint  values,  but  this  effect  can  generally  be  expected  to  be  small. 

The  process  of  making  successive  submissions  has  been  largely  automated 
in  FASTOP-3,  unlike  ASOP-3  where  there  is  a great  deal  of  data  manipulation 
required  in  each  submission.  The  design  produced  by  each  submission  is 
stored  on  a permanent  data  set  and  used  as  input  datp  to  the  next  submission . 
Relatively  little  additional  card  data  are  needed  in  c-scn  sv  .sequent  submission. 

2.3.6  Flexible  Supports 

In  the  resizing  of  a wing  structure  connected  redundantly  to  a fuselage, 
it  is  important  that  the  interactive  forces  at  the  points  of  connection  be 
properly  determined.  A similar  situation  occurs  in  other  cases,  such  as  those 
of  wing  or  tail  surfaces  with  flexible  control  surfaces  or  stores  attached  in 
such  a manner  that  the  interactive  forces  are  statically  indeterminate. 

A capability  of  accommodating  this  type  of  situation  has  been  incorporated 
into  FASTOP-3.  Nodes  that  are  at  points  of  elastic  support  or  points  of 
connection  to  adjoining  elastic  structure  are  considered  to  be  free  to  displace, 
but  additional  stiffnesses  representing  the  constraints  applied  by  the  supports 
or  adjoining  structure  must  be  supplied  by  the  user.  These  additional  stiff- 
nesses are  inserted  at  the  appropriate  locations  in  the  structure  stiffness 
matrix.  In  the  case  of  statically  indeterminate  forces  of  interaction  with 
adjoining  structure,  there  will  generally  be  cross-coupling  between  the  de- 
grees of  f •eedom  at  the  points  of  connection . 


2.4  TRANSFORMATIONS  BETWEEN  MATHEMATICAL  MODELS 

Except  for  some  minor  changes  in  the  required  input  data,  the  transform* 
ation  procedures  in  FASTOP-3  are  identical  to  those  of  FASTOP-1.  Four  dis- 
tinct math  models  are  involved  in  these  transformations,  namely  an  aerodynam- 
ics model,  a weights  model  (for  inertial  loads),  a structures  models,  and  a 
dynamics  model.  Three  transformation  matrices  can  be  created,  which  respec- 
tively transform  loads  from  the  aerodynamics,  weights,  and  dynamics  models 
to  the  structures  model.  The  first  two  matrices  are  needed  for  the  computation 
of  applied  loads,  whereas  the  third  matrix  is  used  in  the  computation  of  a 
dynamics  mcdel  flexibility  matrix.  A comprehensive  treatment  of  the  transforma- 
tion procedures  is  presented  in  Section  S of  Reference  1 . 

2.5  MASS  MATRIX  DEFINITION 

The  procedure  for  obtaining  a dynamics-model  mass  matrix  for  use  in 
vibration  analysis  has  not  changed  in  FaSTOP-3.  There  are  two  alternate 
methods  for  obtaining  this  matrix,  namely  the  "standard"  approach  and  the 
"fully  automated"  approach.  In  the  standard  approach,  the  user  supplies  the 
mass  matrix  for  the  "initial"  design,  i.e..  for  the  design  as  it  exists  prior  to 
any  automated  flutter  resizing.  In  the  fully  automated  approach,  this  matrix 
is  computed  automatically  within  the  program,  based  on  the  weights  of  the 
structural  members  (including  non-optimum  factors,  if  any),  plus  any  fixed- 
mass  items  and  mass-balance  weights  supplied  by  the  user.  Regardless  of 
which  approach  is  used,  the  program  automatically  updates  the  initial  mass 
matrix  to  account  for  the  cumulative  weight  changes  due  to  ali  flutter  and 
strength  resizing  beyond  the  initial  design.  See  Section  6 of  Reference  1 for 
the  details  of  the  procedure. 

2.6  VIBRATION  ANALYSIS 

The  vibration  analysis  procedure  described  in  detail  in  Section  7 of 
Reference  1 applies  to  FASTOP-3  as  well.  Up  to  200  dynamic  degrees  of 
freedom  can  be  accommodated.  Both  cantilever  and  free-free  type  analyses 
can  be  carried  out,  each  using  either  the  "flexibility"  approach  or  the  "sti'f- 
ness"  approach.  For  the  flexibility  approach,  the  user  must  specify  a separate 
dynamics  model,  and  the  eigenvalue  problem  is  then  formulated  in  terms  of  a 


29 


dynamics-model  flexibility  matrix.  The  atiffnesa  approach  is  employed  when 
the  user  does  not  choose  to  provide  a separate  dynamics  model.  In  that 
case,  the  dynamic  degrees  of  freedom  are  taken  to  coincide  with  all  the  struc- 
tures-model  degrees  of  freedom,  and  the  eigenvalue  problem  is  formulated  in 
terms  of  the  structures-model  stiffness  matrix.  The  200- degree-of- freedom 
limitation  mentioned  earlier  obviously  restricts  the  stiffness  approach  to  pro- 
blems in  which  simple  structures  models  are  used. 

2.7  FLUTTER  ANALYSIS 

Except  for  the  improvements  ojscribed  below,  the  flutter  analysis  cap- 
abilities of  FASTOP-3  are  the  same  as  those  of  FASTOP-1,  which  are  described 
in  Section  8 of  Reference  1 . The  program  determines  the  oscillatory  pressures 
and  the  generalized  aerodynamic  forces  for  the  lifting  surface  to  be  analyzed, 
given  a set  of  normal  mode  shapes  and  frequencies.  Generalized  aerodynamic 
forces  are  computed  using  either  the  subsonic  assumed-pressure-function 
procedure  (kernel  function),  the  supersonic  Mach-box  method  or  the  subsonic 
doublet-lattice  procedure.  Then,  given  the  generalized  masses  corresponding 
to  the  modes,  the  program  solves  the  flutter  equation  to  determine  the  lutter 
speed  and  values  of  modal  damping  and  frequency  as  functions  of  air  speed. 
Flutter  solutions  are  obtained  by  use  of  either  the  conventional  k-method  or 
an  improved  version  of  the  p-k  method.  In  the  latter  method,  which  is  used 
for  redesign  purposes,  the  eigenvectors  and  their  associated  row  vectors  are 
determined  at  the  flutter  speed  for  use  in  computing  flutter-velocity  derivatives . 

2.7.1  Improvements  to  the  p-k  Solution  Algorithm 

In  FASTOP-3,  improvements  have  been  introduced  into  the  p-k  solution 
algorithm,  specifically  with  respect  to  root  tracking  and  with  respect  to 
rigid-body  modes.  Modifications  have  been  made  to  the  root  tracking  logic 
because  it  was  discovered  that  the  original  logic,  mder  some  circumstances, 
could  lead  to  incorrect  identification  of  the  critical  root.  Moreover,  with  the 
new  logic,  the  roots  are  generally  determined  after  fewer  iterations,  and 
there  is  more  certainty  that  roots  are  properly  ordered  at  all  velocities. 
Modifications  have  also  been  made  to  permit  inclusion  of  zero-frequency  rigid- 
body  modes  in  the  p-k  flutter  calculations. 


2.7.2 


An  Improvement  to  the  Mach-box  Aerodynamics  Procedure 


Studies  using  the  Mach-box  routines  in  FASTOP-1  showed  that  erratic 
undulations  occurred  frequently  in  the  computed  pressure  distributions. 
Accordingly,  for  FASTOP-3,  these  routines  were  modified  so  as  to  greatly 
improve  the  pressure  distributions  while  retaining  the  essential  flow  discon- 
tinuities. As  a result  of  these  modifications,  some  of  the  text,  equations  and 
figures  in  Subsection  8.2.2  ("Supersonic  Mach-Box  Program")  of  Reference  1 
are  no  longer  applicable.  Accordingly,  a revised  version  of  that  entire  sub- 
section is  presented  on  the  following  pages. 


2.7.2. 1 Revised  Supersonic  Mach-box  Program  for  FASTOP-3 


For  the  supersonic  regime,  the  aerodynamics  program  used  is  a modified 
version  of  the  Mach-box  procedure  described  in  Reference  10.  For  a harmon- 
ically oscillating  planar  surface,  the  pressure  is  related  to  the  velocity  po- 
tential, and  thence  to  the  dowr.wash  distribution,  by 
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where 


x,  i are  streamwise  coordinates 

y,  n are  spanwise  coordinates 

p is  the  differential  pressure  between  the  upper  and  lower  covers  of 

the  surface 

p is  the  air  density 

U is  the  free- stream  velocity 

w is  the  frequency  of  oscillation 

<>  is  the  velocity  potential 

S is  the  lifting-surface  area  bounded  by  the  inverse  Mach  cone 

emanating  from  (x,  y) 

w is  the  complex  downwash  velocity  = U a + iwh  = Uw 

h,  a are  the  deformation  and  slope  of  the  surface 

M is  the  Mach  number 

p = V”2-1 

R = J(x~t)2  - p2  (y-n)2 

” 2 

0^  is  the  differential  pressure  coefficient  = p/(l/2pU  ) 

With  the  exception  of  special  cases,  the  integral  cannot  be  evaluated  in 
closed  form;  hence,  a numerical  approach  is  required.  In  Reference  10,  the 
area,  S,  is  divided  into  elementary  small  rectangular  boxes  having  their  diag- 
onals parallel  to  the  Mach  lines  as  shown  in  Figure  2.6.  The  rectangles  are 
subsequently  converted  to  squares  through  the  coordinate  transformations 

Streamwise  x = x/b, 

l = S/b, 

Spanwise  y = py/b 

n-fin/b,  (2.10) 

where  b is  the  streamwise  dimension  of  a box  and  b/p  is  the  spanwise  dimension. 
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Assuming  the  down  wash  is  constant  over  each  of  these  "Mach  boxes,"  Equation 
(2.9)  can  be  rewritten  as: 
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where 


S.  is  the  j box  within  the  inverse  Mach  cone  emanating  from  (x,  y) 

' ■ f 

s • ■ (S)2  k 

R * ^ (x  - l)2  - (y  - n)2 
b is  the  streamwise  box  size 

Wj  is  the  downwash  on  the  jlh  box 

C.  is  the  jth  pressure  influence  coefficient  for  point  (x,  y),  i.e., 

1 th 

the  pressure  at  the  point  due  to  a unit  downwash  on  the  j box. 

Two  methods  have  been  classically  used  to  perform  the  complex  integration 

over  the  box  area.  The  first,  developed  in  Reference  11,  uses  a mean  value  of 

the  exponential  and  cosine  terms  when  a box  is  far  removed  from  the  point 
~ 2 

(x,  y)  and  a series  expansion  up  to  k when  a box  is  close.  Although  these 
approximations  simplify  the  integration,  they  introduce  significant  errors  when 
the  reduced  frequency,  k.  is  high.  For  a more  exact  evaluation,  a second 
method,  a Bessel-function  series  representation  of  the  integral,  is  presented 
in  Reference  10. 


The  method  used  in  the  current  program  is  different  from  both  of  these 
two  in  that  Gaussian  quadrature  is  used  to  evaluate  the  integral.  With  this 
technique,  singularities  that  occur  in  the  integrand  when  the  box  area  is  cut 
or  touched  by  the  inverse  Mach  cone  emanating  from  the  point  (3,  jf)  - see 
Figure  2.7a  • can  be  accurately  accounted  for  as  described  below.  In  fact,  if 
the  integrand  could  be  represented  exactly  by  the  product  of  a simple 
singular  function,  such  as  1/^/l-y , and  a polynomial  of  order  2N-1,  then 
using  N integration  points  would  produce  the  integral  with  no  error.  Although 
the  true  integrand  cannot  be  presented  exactly  in  the  above  manner,  by 
taking  enough  integration  points  the  error  may  be  made  as  small  as  desired . j 

t 

First,  the  pressure  influence  coefficient  of  Equation  (2.11)  is  rewritten 
by  performing  a change  of  variables: 
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• AT  1 CORNER  • UNIT  SOUARE 

(a)  Type  of  Singularity  and  Integration  Area 
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(b)  Local  Mach  Box  Coordinates  for  Arbitrary  Receiving  Point 


Figure  2.7  Singularities  and  Coordinate  System  for  Mach  Box  Formulation 
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where 


Ax  * x - \ 

Ay  * y - n 

i,  ro  a; the  components  of  tne  distance  between  the  j1*1  box  center 
*nd  the  point  (x.  y)  - see  Figure  2.7b. 

\ « fit  * i?  - (y  - n)2 

ft  s = fi  - «»)2  - (y  - n)  2 

Next,  the  two  single  and  one  double  integrations  are  performed  using  various 
quadrature  formulae  of  Reference  12.  Referring  to  Figure  2.7a,  five  cases  can 
arise  depending  on  how  the  sending  box  is  cut  by  the  inverse  Mach  cone  from 
(x,  y): 

(a)  Box  not  touched:  t i 2,  m $ t - 2. 

(b)  Box  at  apex  of  the  cone:  t - 0,  m = 0. 

(c)  Box  split  by  the  cone:  t - m,  m > 0. 

(d)  Box  touched  at  two  corners  by  the  cone:  i = 1,  m = 0. 

(e)  Box  touched  at  one  comer  by  the  cone:  £>l,m  = *-  l. 

In  case  (a),  since  there  is  no  singularity,  the  following  quadrature 
formula  (Equation  25.4.30,  Reference  12)  may  be  used  for  both  the  spanwise  and 
chordwise  integrations: 


b n 


/ 

f(y)  dy  « « t(yi). 

JB. 

l»l 

where 

b - a 

b ♦ a 

yi  * 2 *i 

♦ 2 

Xj  is  the  i1*1 

zero  of  the  Legendre  polynomial,  P] 

H.,  the  weighting  function  at  the  i1^  quadrature  point, 
= 2/(1  - x*  ) ^ (Xj  )J  2 


(2.14) 
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In  case  (b),  the  integration  is  to  be  performed  only  on  the  triangular 
area  of  a quarter-box.  Hence,  reversing  the  order  of  integration  in  Equation 

V Ax 

(2.13),  the  limits  of  the  double  integral  become  / / , while  the  limits 

A •'O  •'-Ax 

of  the  single  integral  are  / . Since  R = 0 at  the  limits  Ax  and  -Ax,  and 


since  Rj^  * ^-1/2  = 0 81  the  limits  1/2  and  -1/2,  the  spanwise  integrals  have 
singularities,  and  a different  quadrature  formula  (Equation  25.4.39,  Reference 
12)  is  used: 


(2.15) 


where 


y = L.:  -«  x + k.  JL« 

yi  2 * 2 

Xj  = cos  (2i  - 1) 

H.  = n/n. 

After  the  spanwise  integration  is  performed,  Formula  (2.14)  is  used  in  the 
chordwise  direction. 

For  case  (c),  the  integration  is  done  over  a triangular  half-box  where 

-<♦<»  -Ax 

the  limits  of  the  integrals  become  / I and  I Since  the  span- 

' J 1-* 

wise  integrals  here  have  singularities  at  the  upper  limits  only,  the  appropriate 
quadrature  formula  (Equation  25.4.37,  Reference  12)  is: 


- - (2n) 

Vb  • v./  Hi 

1-1 


(2.16) 


where 


y.  « a*  (b-a)  (1-x2.) 

Xj  is  the  i1*1  positive  2ero  of  the  Legendre  polynomial,  P^. 

H^2n^  are  the  Gaussian  weights  of  order  2n. 

Again,  Formula  (2.14)  is  used  in  the  chord  wise  integration. 

For  case  (d),  an  integration  is  first  performed  over  the  triangular 
shaded  area  in  Figure  2.7a  (comprising  boxes  b,  c,  and  d)  using  the  quad- 
rature formulae  for  case  (b).  Then  by  subtracting  the  previously  derived 
pressure  influence  coefficients  for  boxes  b and  c,  the  desired  integral  over  d 
is  achieved. 

Finally,  for  case  (e),  the  integral  for  an  aggregate  area  (see  dotted  area 
in  Figure  2.7a)  consisting  of  a triangle  on  the  Mach  line  and  the  subject  area 
is  determined;  and  from  it  the  integral  for  a triangle  is  subtracted.  When 
the  spanwise  integrations  are  performed,  Equation  (2.16)  is  employed,  while 
Equation  (2.14)  is  used  in  the  chordwise  direction. 

For  case  (a)  and  all  chordwise  integrations,  the  present  program  uses 
six  integration  points  for  the  quadrature.  In  the  spanwise  integrations,  six 
points  are  used  when  k < \ *£)2,  while  twelve  points  are  employed  when 

k>M>2 

K - 3 *m;  • 

At  a given  Mach  number  and  reduced  frequency,  the  pressure  influence 
coefficients  are  functions  of  only  £ and  m - the  separation  between  the  sending 
and  receiving  box  centers.  Consequently,  influence  coefficients  are  computed 
by  the  above  formulae  only  once  for  each  sdr.issible  £,  m pair  (£>0,m<£) 
and  are  used  repeatedly  where  needed. 

The  pressure  on  any  box  is  a function  of  the  downwash  of  only  those 
boxes  within  the  inverse  Mach  cone  emanating  from  its  box  center.  For  a 
surface,  the  edges  of  which  are  all  supersonic,  the  pressure  is,  furthermore, 
only  influenced  by  boxes  on  the  planform.  If  any  of  the  surface  edges  are 
subsonic,  however,  there  are  regions  adjacent  to  these  edges  which  do  affect 
the  pressure  of  some  areas  on  the  planform.  To  account  for  this  effect,  the 
concept  (Reference  13)  of  a permeable  "diaphragm"  is  introduced  in  these 
regions.  This  permeable  sheet  does  not  alter  the  flow  nor  can  it  sustain 
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pressure.  It  is  bounded  by  the  surface  edge  and  the  Mach  lines  emanating 
from  the  comers  of  the  lifting  surface  - see  Figure  2.8. 

The  relationship  between  the  pressure  and  downwash  on  the  combination  of 
the  lifting  surface  and  the  diaphragm  area  can  be  written: 


pw  I _ fcwwl  c wdI  (Ww| 

’pd(  [j5 Dw!  ^ DDJ  1 1 D j 


(2.17) 


where 


rW 


'WW 


'WD 


'DW 


'DD 


w 


W 


wr 


is  the  pressure  on  the  wing  boxes 
is  the  pressure  on  diaphragm  boxes 

are  the  influence  coefficients  giving  pressures  on  the  wing  boxes 
due  to  downwash  on  the  wing  boxes 

are  the  influence  coefficients  giving  pressures  on  the  wing  boxes 
due  to  the  downwash  on  the  diaphragm  boxes 

are  the  influence  coefficients  giving  pressures  on  the  diaphragm 
boxes  due  to  downwash  on  the  wing  boxes 

are  influence  coefficients  giving  pressures  on  the  diaphragm  boxes 
due  to  downwash  on  the  diaphragm  boxes 

is  the  known  downwash  on  the  wing  boxes 

is  the  unknown  downwash  on  the  diaphragm  boxes. 


Since  the  pressure  on  any  diaphragm  box  is  zero,  then 


[CDw]  {Ww}*  [CDd]  {Wd}=  {°}’ 


and  the  unknown  diaphragm  downwash  can  be  evaluated  by 

,-l 


{wd}  = - £cdd]  [cdw]  {Ww} 


(2.18) 


(2.19) 
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Substituting  this  result  into  Equation  (2.17)  yields  the  final  expression  for 
the  pressures  on  the  wing  boxes: 


{pw}  * M M 

-i 

A1C1  s M-  [cwd][cdq1  [cdw] 


(2.20) 

(2.21) 


For  maximum  computer  efficiency,  the  actual  calculation  of  the  pressures. 
Pw,  is  generally  performed  in  a different  manner  from  that  implied  by  Equation 
(2.20).  The  calculation  of  the  aerodynamic  influence  coefficient  matrix,  (A1C), 
leads  to  either  extenv.ve  use  of  core  in  storing  matrices  or  a large  number  of 
I/O  operations  if  the  .matrices  are  stored  on  data  devices.  If  there  is  no  need 
for  saving  the  {A1C)  array  for  future  use,  the  machine  operations  can  be  appre- 
ciably reduced  by  computing  the  pressures  as  follows: 


(1) 

\y\  - 

% 

i 

D 

O 

(2) 

M = 

(CDdI!  J X 1 
IcwdI  I'M 

(2.22) 

(3) 

{ z i = 

(4) 

1 pw!  ’ 

lcww]  l"wl  • lZl 

In  this  way.  core  storage  requirements  are  minimized,  since  only  a vector  need 
be  stored  in  going  fron  one  step  to  another. 

Once  the  oscillatory  pressures  are  computed,  the  normalized  generalized 
aerodynamic  forces  are  computed  by 


I - ff  y'  -X 

*s  JJ..,  ’ ** 


(2.23) 


Since  the  boxes  are  assumed  to  be  very  small,  the  integration  can  be  replaced 
by  a summation  over  every  box  on  the  planform: 


’ 23f  4-  r 


(2  24) 


A 


) 


is  the  deflection  of  the  j1*1  box  in  the  r1*1  mode 
is  the  pressure  coefficient  of  the  j1^  box  in  the  s1*1  mode 
is  the  srea  of  the  j1*1  surface  box . 

2 

It  should  be  noted  that  each  box  area  is  b /p.  The  program  automatically 
establishes  the  gridwork  of  boxes:  From  a user  specified  number  of  boxes 

desired,  the  program  calculates  the  box  si2e  necessary  for  the  boxes  to  cover 
the  planform  and  diaphragm  and  to  align  with  the  inboard  and  outboard 
planfortn  edges.  Consequently,  no  boxes  overhang  the  planform  side-edges. 
(See  Figure  2.8.)  Each  box  is  designated  as  either  a wing  box  or  diaphragm 
box  depending  on  whether  its  center  lies  on  or  off  the  wing,  respectively. 
In  general,  boxes  do  overhang  the  leading  and  trailing  planform  edges 
causing,  in  effect,  a jagged  representation  of  these  edges.  Although  for 
most  configurations,  this  jaggedness  has  been  found  to  have  little  effect  on 
the  accuracy  of  the  generalized  aerodynamic  forces  for  simple  mode  shapes, 
providing  that  the  box  grid  is  not  too  coarse,  it  does  produce  erratic  undu- 
lations in  the  computed  pressures . 

A recent  analytical  refinement  (see  Reference  14)  remedies  this  short- 
coming of  the  basic  method,  and  produces  greatly  unproved  supersonic  oscilla- 
tory pressure  distributions  with  only  a nonimal  increase  in  computer  time. 
This  refinement,  which  has  been  incorporated  into  FASTOP-3.  is  described  as 
follows : 

Any  box  that  is  crossed  by  a planform  edge  is  considered  to  consist  of 
two  identical  overlapping  boxes,  one  being  treated  as  if  it  were  on  the  plan- 
form  and  the  other  as  if  it  were  on  the  diaphragm.  Let  the  pressures,  down- 
washes,  and  influence  coefficients  for  such  boxes  be  denoted  by  the  subscripts 
WS  (wing,  shared)  and  DS.  (diaphragm,  shared).  Now,  since  the  area  of  the 
planform  encaptured  by  a box,  WS{,  is  only  a fraction,  A^g  /A,  of  the  full  box 

area  (see  box  A,  Figure  2.8),  the  contribution  to  the  pressure  at  any  other 
box  due  to  a signal  originating  from  this  box  should  be  scaled  by  this  ratio. 
Although  similar  reasoning  might  be  applied  to  the  box  DS^,  it  has  been  demon- 
strated in  sample  calculations  that  this  scaling  should  only  be  applied  to  the 
wing-shared  boxes,  incorporating  these  ideas.  Equation  (2.17)  is  replaced  by 
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where  the  factors  of  1/A  ar^  omitted  since  the  area  df  a full  box  in  the 
transformed  Mach-box  system  is  identically  unity.  The  2ero  submatrices  in 
this  equation  are  necessary  to  assure  non-zero  leading-edge  pressures.  For 
further  explanation,  refer  to  Reference  14. 


An  additional  distinction  in  handling  boxes  is  required  to  properly  treat 
trailing  edges.  For  a typical  box  on  a supersonic  trailing  edge,  the  inverse 
Mach  cone  fi-om  a point  at  the  center  of  the  box  encompasses  the  same  wing 
area'.  A/4,  as  does  any  upstream  whole  wing  box  (see  shaded  area  in  box  B. 
Figure  2.8).  Since  the  downwash  on  such  a box  affects  the  pressure  on  almost 
no  other  box.  the  area  ratio  js  not  applied  to  the  influence  coefficient  of 
this  box  For  a box  on  a subsonic  trailing  edge,  the  wing  area  encompassed  by 
the  inverse  Mach  cone  is  still  generally  a much  larger  percentage  of  the  total 
Mach  cone  area  than  is  indicated  by  the  ratio.  A^gi  hence,  the  influence 
coefficient  for  this  box  is  also  not  scaled. 


Results  obtained  by  the  above  formulation  require  some  interpretation. 
Because  of  the  jagged  representation  of  the  leading  edge,  the  centers  of  some 
boxes,  such  as  box  C.  Figure  2.8.  lie  physically  ahead  of  the  true  planform; 
yet.  non-zero  pressures  arc  computed  for  such  boxes.  This  is  correct,  since 
the  calculated  pressure  coefficient.  f*Wg.  is  applied  over  the  entire  box  area 
and.  certainly,  the  pressure  over  A^g  is  non-zero.  Now.  if  the  force  Pwg  • A 
is  applied  to  only  the  correct  area.  A^g,  the  true  pressure  coefficient  for 
that  area  would  be  Pws  . A/A^g.  or,  since  A = 1 , P^S^^WS'  However,  since  the 
contribution  from  such  shared  boxes  to  the  generalized  forr*  is  given  by 


^ " (fAws]  * ,Pws/Aws’)  * lhws* 

*PWS*  ' 


(2.26) 


the  pressures  computed  by  solving  Equation  (2.25)  rather  than  the  weighted 
values.  P,g/A^g.  are  actually  used  in  the  computer  program. 
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Results  of  numerous  examples  in  Reference  14  show  that  the  pressure 
distributions  obtained  from  the  present  method  are  substantially  better  than 
those  computed  by  the  basic  method. 

When  a highly  swept  surface  is  analyzed  for  a relatively  low  Mach  number, 
the  number  of  forward  diaphragm  boxes  can  become  so  great  as  to  cause  an 
appreciable  increase  in  computing  time.  However,  the  down  wash  in  the  dia- 
phragm region  decreases  very  ...pidly  in  the  forward  streamwise  direction. 
To  save  computing  time,  the  present  program  takes  this  rapid  decay  of  dia- 
phragm downwash  into  account  and  allows  the  user  to  request  a box- 
elimination  option  whereby  the  diaphragm  boxes  are  ignored  forward  of  a 
user-specified  distance  ahead  of  the  leading  edge. 

In  the  program,  provision  is  made  for  computing  aerodynamic  force 
coefficients  and  center-of-pressure  locations.  The  user  may  use  this  facility 
to  compare  known  steady-state  data  with  computed  values  to  determine  the 
number  of  boxes  required  for  a satisfactory  solution.  Another  approach  is  to 
vary  the  number  of  boxes  and  look  for  convergence  in  the  stability  coefficients. 
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2.8  RESIZING  FOR  COMBINED  FLUTTER  AND  STRENGTH  REQUIREMENTS 

Except  for  a few  details  regarding  the  treatment  of  composite  members 
the  procedure  described  in  Section  9 of  Reference  1 for  effecting  interactive 
strength/flutter  redesign  still  applies  in  FASTOP-3.  However,  as  this  pro- 
cedure is  a major  feature  of  the  FASTOP  system,  it  is  set  forth  again  below. 
The  reader  should  consult  the  original  reference  for  the  details  of  the  re- 
quired computational  operations,  e.g.  , computation  of  flutter-velocity  derivatives. 

2.8.1  Interactive  Strength/Flutter  Resizing  Procedure 

Starting  with  the  use  of  the  Structural  Optimization  Program  (SOP),  the 
structure  is  sized  to  satisfy  its  strength  requirements  with  a fully  stressed 
design.  Two  categories  of  elements  exist  after  this  first  step;  specifically, 
each  element  either  is  fully  stressed  (i.e.,  "strength-critical")  or  is  at  its 
prespecified  minimum  gage  (as  dictated,  for  example,  by  manufacturing  con- 
siderations). The  next  step  uses  the  Flutter  Optimization  Program  (FOP)  to 
resize  structural  elements  and/or  mass-balance  design  variables  to  increase 
the  surface's  critical  flutter  speed.  None  of  the  structual  elements  are 
permitted  to  be  reduced  in  size  in  this  initial  FO?  step  since,  upon  entering 
the  step,  all  of  these  elements  were  already  either  strength-critical*  or  at 
prespecified  minimum  gage.  Those  structural  elements  that  are  increased  in 
size  in  this  step,  plus  any  mass-balance  variables  present  in  the  design,  con- 
stitute a third  category  of  elements,  namely,  "flutter-critical"  elements,  i.e., 
elements  whose  gages  are  dictated  by  flutter-speed  requirements. 

The  resizing  of  some  structural  elements  during  the  first  FOP  step  may 
cause  a significant  redistribution  of  internal  loads  within  the  structure, 
thereby  modifying  gage  requirements  for  strength  considerations.  Accordingly, 
in  the  next  step,  SOP  is  re-entered  for  the  purpose  of  a "strength  update". 
This  is  the  first  attempt  to  achieve  a minimum  weight  design  that  accounts  for 
strength/flutter  interaction.  In  this  second  SOP  step,  the  flutter-critical 
elements  (and,  of  course,  the  minimum-gage  elements)  are  not  permitted  to  be 


*In  FASTOP- 1,  member  stress  ratios  were  passed  from  SOP  to  FOP  where  they 
served  to  define  minimum  size  requirements  for  strength.  In  FASTOP-3,  these 
minimums  are  passed  directly,  i.e.,  stress  ratios  are  no  longer  used  in  FOP. 
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resized  downward,  but  the  strength-critical  elements  are  free  to  be  sized  up- 
ward or  downward.  After  resizing,  the  various  elements  are  reclassified  into 
strength-critical,  flutter-critical,  and  minimum-gage  categories.  Elements 
may,  of  course,  change  categories.  For  example,  an  element  would  shift  from 
the  flutter-critical  to  the  strength-critical  set  if  it  had  been  resized  upward 
to  satisfy  the  modified  strength  requirements.  Likewise,  an  element  which 
was  previously  strength-critical  might  be  resized  downward  to  the  point  where 
it  enters  the  minimum-gage  category . 

At  this  stage,  FOP  is  entered  for  the  second  time  and  the  interactive 
strength/flutter  redesign  process  continues.  As  in  the  first  FOP  step,  ele- 
ments  in  the  strength-critical  and  minimum-gage  categories  cannot  be  sized 
downward.  On  the  other  hand,  there  now  exists  a set  of  elements,  namely 
the  flutter-critical  elements,  which  are  free  to  be  resized  in  either  direction; 
if  such  an  element  is  sized  downward,  however,  its  gage  is  not  permitted  to 
f-’ll  below  the  values  required  by  the  last  strength  analysis  or  by  minimum- 
gage  considerations.  After  resizing,  the  elements  again  are  reclassified  into 
the  three  basic  categories. 

Subsequent  interactive  application  of  the  two  programs  proceeds  in  a 
manner  similar  to  the  second  SOP  and  FOP  steps  until  the  process  is  suffi- 
ciently converged.  The  final  design  will  consist  of  a set  of  flutter-critical 
elements  which  have  nearly  uniform  flutter- velocity  derivatives,  a set  of 
strength-critical  elements  which  are  fully  stressed,  and  a set  of  elements 
which  are  at  prespecified  minimum  gages . 

2.8.2  Flutter  Resizing  of  Composite  Members 

Flutter  resizing  of  composite  members  involves  a few  considerations  that 
have  no  counterparts  for  metallic  members,  namely  the  "rounding"  and 
"balancing"  requirements.  These  considerations  are  discussed  below. 

2 . 8 . 2 . 1 Rounding  Requirement 

For  a composite  member,  the  design  variables  are  the  number  of  laminae 
in  each  of  the  layers.  When  the  flutter  resizing  equation  is  applied  to  such 
a variable,  the  new  value  of  the  variable  will  generally  nol  be  an  integral 
number  of  laminae.  Accordingly,  in  FASTOP-3,  it  has  been  decided  to  round 
the  new  value  to  the  nearest  integer. 
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2. 8. 2. 2 Balanced-Design  Requirement 

In  the  conventional  use  of  FASTOP-3  for  strength/flutter  redesign, 
either  balanced-laminate  resizing  or  unbalanced-laminate  resizing  is  used 
consistently  throughout  the  entire  redesign  process,  i.e.,  for  both  strength 
and  flutter  resizing.  The  user  specifies,  by  means  of  clue  data  in  the  first 
SOP  pass,  which  of  the  two  possible  types  of  resizing  is  to  be  performed  for 
all  subsequent  strength  and  flutter  redesign.  There  is  a situation,  however, 
in  which  departure  from  this  consistent  resizing  approach  may  be  desirable. 
Consider  the  case  in  which  a balanced-laminate  fully  stressed  design  has  been 
generated  in  an  initial  SOP  pass.  If  this  fully  stressed  design  is  subsequently 
found  to  have  a marginal  flutter-speed  deficiency,  the  user  might  wish  to 
resolve  the  deficiency  by  merely  adding  a "patch"  of  unbalanced  material  to 
the  design.  Accordingly,  FOP  contains  a data  clue  to  make  this  possible.  If 
that  clue  is  exercised,  however,  the  user  logically  should  not  expect  to  call 
for  further  strength  resizing  in  any  subsequent  pass  through  SOP,  since  that 
program  would  again  generate  a balanced  design;  that  is,  SOP  should  thereafter 
be  used  to  analyze  designs  only. 

When  called  for,  balanced-laminate  flutter  resizing  is  accomplished  by  (a) 
computing  the  true  flutter-velocity  derivatives  for  the  layers  to  be  balanced, 
(b)  averaging  these  derivatives,  and  (c)  using  the  average-derivative  value 
when  applying  the  resizing  equation  to  each  of  the  pertinent  layers.  This 
procedure  is  also  used  for  balanced- laminate  deflection  resizing  in  SOP. 

2.9  USE  OF  FASTOP  FOR  INTEGRATED  ANALYSIS  AND  DESIGN 

Considerations  on  the  use  of  FASTOP  for  interactive  strength/flutter 
redesign,  as  described  in  Section  10  of  Reference  1,  also  apply  to  FASTOP-3. 
The  user-oriented  material  presented  in  that  reference  is  intended  to  guide 
the  user  through  the  entire  strength/flutter  redesign  process.  To  aid  the 
user,  numerous  suggestions  are  provided  based  on  experience  obtained  during 
the  performance  of  the  demonstration  problems. 


Section  3 


EXAMPLES  OF  PROGRAM  APPLICATIONS 


3.1  INTRODUCTION 

FASTOP-3  has  been  applied  to  two  different  structures  to  demonstrate 
the  program's  capability  to  obtain  near-optimun  designs  for  metallic  and  com- 
posite lifting  surfaces  subjected  to  both  strength  and  flutter-speed  require- 
ments. The  two  cases  selected,  namely  the  Intermediate-Complexity  Wing  and 
the  All-Movable  Stabilizer,  represent  levels  of  structural  detail  that  would 
exist  at  different  stages  of  the  design  process.  Specifically,  the  coarse  model 
of  the  wing  is  representative  of  the  preliminary  design  phase,  whereas  the 
more  detailed  stabilizer  model  is  typical  of  what  might  be  used  during  the 
later  stages  of  design.  The  two  models  and  the  results  of  the  interactive 
strength/flutter  redesign  studies  are  described  in  the  following  subsections. 
In  addition,  results  are  presented  for  resizing  of  the  Intermediate-Complexity 
Wing  to  satisfy  deflection  constraints. 

3.2  INTERMEDIATE-COMPLEXITY  WING 
3.2.1  Mathematical  Models 

Figure  3.1  shows  the  aerodynamic  planform  and  primary  structural  arrange- 
ment of  the  Intermediate-Complexity  Wing.  The  structures  model  consists  of  a 
cantilevered,  symmetric,  two-cell  box  beam,  having  aluminum  substructure 
and  graphite/epoxy  composite  cover  skins  in  a 0°/±45°/90°  layup.  The  zero- 
degree  fibers  are  aligned  with  the  center  spar  of  the  box.  Membrane  elements 
are  used  to  model  the  skins;  shear  panels  represent  the  rib  and  spar  webs; 
and  bar  elements  are  introduced  between  upper-  and  lower-cover  node  points. 
The  model  has  234  degrees  of  freedom  and  158  members,  64  of  which  are 
composite  members  requiring  4 design  variables  each . 

The  dynamics  model  of  the  Intermediate-Complexity  Wing  is  shown  in 
Figure  3.2.  There  are  39  dynamics  nodes,  each  of  which  lies  in  the  wing 
mid-plane  directly  between  upper  and  lower  cover  structural  node  points. 
Every  dynamics  node  is  permitted  to  undergo  out-of-plane  displacement.  In 
addition,  nodes  along  the  perimeter  of  the  wing  also  have  rotational  degrees 
of  freedom  about  either  the  front  spar,  rear  spar,  tip  rib,  or  pitch -axis 
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directions.  These  rotational  freedoms  are  needed  to  properly  account  for  the 
mass  and  inertia  of  the  overhung  structure.  In  all.  the  dynamics  model  has 
58  degrees  of  freedom. 

3.2.2  Results  of  Strength/Flutter  Redesign  Study 

Using  two  applied  load  conditions  (one  subsonic  and  one  supersonic)  the 
Structural  Optimization  Program  (SOP)  was  employed  to  obtain  fully  stressed 
designs  (5  FSD  cycles)  for  both  the  unbalanced*  and  balanced-laminate  options; 
balancing  was  with  respect  to  the  ±45°  layers.  The  resulting  unbalanced 
design  weighed  37.7  lb,  14.3  lb  of  which  went  into  each  of  the  two  identical 
covers.  For  the  balanced  design,  the  weight  was  41.6  lb,  each  cover  weighing 
16.3  lb. 

The  Flutter  Optimization  Program  (FOP)  was  then  used  to  perform  vibra- 
tion and  flutter  analyses  for  both  the  unbalanced  and  balanced  fully  stressed 
designs.  The  flutter  analyses  utilized  the  doublet- lattice  procedure  for  Mach 
0.8,  sea  level,  with  the  six  lowest-frequency  normal  modes  of  vibration.  For 
the  unbalanced  design,  the  flutter  speed  was  712  KEAS;  the  balanced  design 
had  a flutter  speed  of  780  KEAS. 

Two  parallel  interactive  strength/flutter  redesign  studies  were  then  ini- 
tiated to  raise  the  flutter  speeds  of  both  designs  to  925  KEAS  (1.3  x 712). 
Flutter  resizing  was  restricted  to  the  composite  cover  skins  only,  but  all 
elements  were  candidates  for  strength  resizing.  Figure  3.1  shows  a history 
of  the  resizing  steps  for  both  balanced  and  unbalanced  laminates.  In  each 
case,  an  efficient  final  design  was  obtained  in  6 combined  strength/flutter 
resizing  steps.  The  first  3 steps  in  each  study  were  used  to  raise  the  flutter 
speed  from  the  FSD  value  to  the  vicinity  of  the  specified  value  of  925  KEAS. 
Additional  steps  were  then  performed  to  approach  a converged  design. 

As  shown  in  Figure  3.3,  the  unbalanced-laminate  final  design  had  a 
structural  weight  of  approximately  44.0  lb,  or  about  6.3  lb  more  than  the 
starting  fully  stressed  design.  The  detailed  results  indicated  that  this  weight 
increment  was  associated  almost  entirely  with  the  cover  skins,  and  that  approx- 
imately 7.5  lb  had  been  added  to  achieve  the  increased  flutter  speed  while 
about  1.2  lb  had  been  removed  due  to  the  decreased  stress  levels  resulting 
ttx>M  strength/flutter  interaction.  The  unbalanced- laminate  layups  for  the 
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fully  stressed  design  and  for  the  final  design  are  shown  in  Figure  3.4.  Note 
the  large  increase  in  the  thickness  of  the  0°  layers  in  the  trailing-edge 
panels  near  the  root.  The  weight  increase  in  this  zone  accounted  for  almost 
70%  of  the  total  weight  increment . 

The  results  of  the  balanced-laminate  redesign  were  qualitatively  similar 
to  the  unbalanced  case  in  that  the  weight  increment  was  again  primarily  in 
the  cover  skins  and  involved  a moderate  amount  of  strength/flutter  interaction. 
Figure  3.5  shows  the  balanced- laminate  layups  for  the  fully  stressed  design 
and  the  final  design. 

3.2.3  Results  of  Strength/Deflection  Redesign  Study 

The  deflection-constraint  resizing  study  originally  performed  on  the 
Intermediate-Complexity  Wing  to  demonstrate  the  capability  of  the  ASOP-3 
program  (Reference  4)  has  been  repeated  here  to  verify  the  same  capability 
in  FASTOP-3. 

The  fully  stressed  design  for  the  unbalanced  layup,  discussed  in  the 
preceding  subsection,  was  examined  from  the  point  of  view  of  streamwise-twist 
distribution  along  the  wing's  span  for  the  subsonic  loading  condition  - the 
twist  angle  being  based  simply  on  the  difference  in  vertical  displacements 
between  the  forward  and  aft  wing  spars  along  a streamwise  chord.  This  twist 
distribution  is  shown  by  the  upper  curve  in  Figure  3.6.  It  is  interesting  to 
note  that  the  forward  center  of  pressure  of  the  subsonic  loading  distribution 
causes  sufficient  nose-up  twisting  to  overpower  the  usual  nose-down  twisting 
(washout)  that  generally  occurs  in  swept  metallic  wings. 

To  illustrate  a potential  application  of  the  deflection-constraint  resizing 
capability  of  FASTOP-3,  it  was  uecided  to  attempt  to  "tailor"  the  design  to 
achieve  a prescribed  streamwise-twist  distribution  for  the  subsonic  loading 
condition  that  would  offer  improved  aerodynamic  performance  through  increased 
lift-to-drag  ratio.  Twist  angles  at  two  wing  stations  were  then  established  as 
targets,  these  being  -2.0°  (washout)  at  a selected  inboard  station  and  -2.5° 
at  the  most-outboard  rib  station  (see  Figures  3.1  and  3.6). 

Resizing  in  the  deflection-constraint  mode  was  accomplished  in  two  stages. 
The  approach  was  to  divide  the  structure  in  two  regions,  as  indicated  by  the 
bold  separating  line  in  Figure  3.1.  In  the  first  resizing  stage,  only  the 
composite  cover  skin  elements  in  the  inboard  region  were  permitted  to  be 
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issued  in  the  deflection  part  of  a resizing  cycle,  to  meet  the  inboard  station 
twist-angle  requirement  In  the  second  stage,  only  out  board- region  cover 
skin  elements  Mere  allowed  to  be  resized,  to  achieve  the  desired  outboard- 
station  l vast  angle.  In  ooth  stages,  however,  all  elements  were  permitted  to 
be  resized  if  they  were  strength  critical.  This  two-stage  approach  was  based 
on  the  concept  that,  for  high-aspect-ratio  cantilevered  surface the  resizing 
of  elements  outboard  of  a particular  station  should  have  little  influence  on  the 
deflections  at  that  station. 

The  first  stage  of  resizing  in  the  deflection -constraint  mode  started  with 
the  fully  stressed  design.  Convergence  to  the  desired  twist  angle  at  the 
inboard  station  was  achieved  in  eight  : eps,  with  the  overall  resulting  twist 
distribution  as  shown  by  the  dashed  curve  in  Figure  3.8.  Figure  3.7  sum- 
marizes the  resizing  history  in  this  mode,  in  terms  of  inboard- station  twist 
angle  versus  total  structural  weight,  after  the  strength- resizing  part  of  each 
cycle . 

In  the  second  stage  of  deflection-constraint  resizing  all  starting  gages 
were  taken  as  those  of  the  final  design  in  the  first  stage.  For  elements  in 
the  inboard  region,  these  starting  gages  were  also  treated  as  mmimums.  to 
prevent  removal  of  material  that  was  previously  introduced  to  meet  the 
inboard-station  twist  requirement.  It  should  be  noted  that  the  introduction  of 
these  starting  and  minimum  gages  is  fully  automated  in  FASTOP-3  by  making 
use  of  a stored  member  data  file  generated  in  the  first  resizing  stage.  This 
is  different  from  ASOP-3.  where  considerable  hand  manipulation  of  member  data 
is  necessary 

Convergence  to  the  desired  outboard- station  twist  angle  required  only 
two  cycles,  with  only  a very  small  additional  weight  increase.  The  final  twist 
distribution  after  this  second  deflection-constraint  resizing  is  shown  in  Figure 
3 6.  and  a summary  of  results  for  all  stages  of  resizing  is  presented  in  Table 
3 1.  It  should  be  pointed  out  that  the  small  differences  between  the  target 
and  accepted  twist  angles  are  due  mainly  to  limits  imposed  by  the  practical 
requirement  for  rounding  layers  to  integral  "lumbers  of  laminae  Figure  3 8 
displays  the  wing  cover  layups  for  the  initial  fully  stressed  design  and  the 
final  combined  strength-and  deflection-constrained  design 
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Figura  3.7  Intarmadiata-Complaxlty  Wing  - Ratizlng  Hiltory  from  Fully  Strauad  Datign 
to  Daftaetlon-Conitrainad  Daiign  for  Inboard  Station 
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Fiflura  3.8  Covtr-Skin  Layups  for  IntarmadiataCompiaxity  Wing  (Unbaiancad-Laminata 
Case)  Rasuad  to  Satisfy  Strangth  and  Daflaction  Constraints 
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TABLE  3.1.  SUMMARY  OF  RESULTS  FOR  STRENGTH  AND  DEFLECTION- 
CONSTRAINT  RESIZING  OF  INTERMEDIATE-COMPLEXITY  WING 


Constraint 

Mode 

Cycles  To 
Convergence 

(1) 

6 Inbd 

(2) 

0 Outbd 

Structural 

Weight 

Stress 

5 

♦0.50° 

♦ 1 .09° 

37.7  lb 

First  Deflection 

8 

-2.01 

-1.93 

GO. 6 

Second  Deflection 

2 

-2.03 

-2.51 

61.2 

^Desired  Value  = -2.00°  (Leading  Edge  Down) 
( 21 

' Desired  Value  = -2.50°  (Leading  Edge  Down) 


3.3  ALL-MOVABLE  STABILIZER 
3.3.1  Mathematical  Models 

The  aerodynamic  planform  and  structural  idealization  of  the  all-movable 
stabilizer  is  shown  in  Figure  3.9.  Cover  skins,  modeled  as  membrane  elements, 
are  of  composite  construction,  except  for  small  2ones  of  metallic  construction 
(the  "splice-plates")  in  the  vicinity  of  the  support  points.  A hybrid  laminate 
was  selected  to  demonstrate  the  capability  of  the  program.  This  laminate  is 
constructed  of  boron/epoxy  in  the  0°  (spanwise)  direction,  and  graphite/ 
epoxy  in  the  ±45°  and  90°  directions.  Titanium  is  used  in  the  splice  plate. 
The  substructure,  which  is  aluminum  honeycomb  material,  is  modeled  as 
spanwise  and  ribwise  shear  panels  having  stiffness  equivalent  to  the 
honeycomb.  Posts  are  also  present  between  upper-  and  lower-cover  node 
points.  In  addition,  there  are  metallic  redistribution  ribs  and  a few 
metallic  spar  webs  in  the  vicinity  of  the  splice  plates . The  details  of  the 
pivot  restraints  and  the  actuator  arm  are  omitted  from  the  figure  to  preserve 
clarity  of  presentation.  It  should  be  noted,  however,  that  the  inner  and 
outer  stabili2er-to-pivot  supports  are  modeled  in  the  mid-plane  of  the  surface 
at  points  respectively  enclosed  by  node-sets  481,  483,  463,  461  and  379,  381, 
349,  347.  In  all,  the  model  has  1172  degrees  of  freedom  and  891  finite 
elements.  Each  cover  has  162  members,  all  but  15  of  which  are  composite 
elements  requiring  4 design  variables  each. 
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Figure  3.10  shows  the  dynamics  model  of  the  all-movable  stabilizer.  F.ach 
of  the  73  dynamics  nodes  lies  in  the  mid-plane  of  the  structure  directly  be- 
tween upper-  and  lower-cover  structural  node  points.  Out-of-plane  displacement 
is  permitted  at  every  dynamics  node.  In  addition,  to  account  for  the  effects 
of  the  overhung  structure,  19  nodes  along  the  front  and  rear  spars  have  ro- 
tational degrees  of  freedom  about  either  the  front  or  rear  spar  directions. 
The  model  therefore  has  a total  of  92  degrees  of  freedom . 

3.3.2  Results  of  Strength/Flutter  Redesign  Study 

In  the  redesign  study  of  the  all-movable  stabilizer,  resizing  to  satisfy 
both  strength  and  flutter  requirements  was  restricted  to  the  elements  in  the 
cover  skins  only.  All  other  elements  in  the  structures  model  were  effectively 
excluded  from  the  entire  redesign  process  because,  for  each  of  those  elements, 
the  minimum  and  maximum  gages  were  set  equal  to  the  starting  gage.  As  a 
result,  the  design  problem  involved  1206  active  design  variables,  all  but  30  of 
which  were  associated  with  the  layers  of  the  composite  elements  in  the  skins. 
T^e  remaining  30  variables  were  associated  with  the  metallic  splice-plate 
elements . 

For  the  stabilizer  structures  model  described  previously,  fully  stressed 
designs  were  obtained  in  5 cycles  for  both  the  balanced-  and  unbalanced- 
laminate  options;  balancing  was  with  respect  to  the  ±45°  layers.  Three  applied 
loading  conditions  were  used,  two  of  which  (M  = 0.8,  S.L.  and  M = 1.3, 
10,000  ft.)  were  generated  within  the  automated  load  analysis  module,  ALAM. 
The  unbalanced  ^nd  balanced  designs  weighed  207  lb  and  214  lb,  respectively. 
Figure  3.11  shows  the  upper-cover  layup  of  the  unbalanced- laminate  case. 
Note  that  all  layers  contain  even  numbers  of  laminae  to  enforce  laminate 
designs  having  mid-plane  symmetry.  This  feature  was  obtained  by  specifying 
double  thicknesses  for  the  laminae  in  the  materials  input  data.  Lower-  and 
upper-cover  layups  were  almost  identical. 

The  Flutter  Optimization  Program  (FOP)  was  then  used  to  perform  vibration 
analyses  and  Mach-box  flutter  analyses  for  both  the  unbalanced-  and  balanced- 
laminate  fully  stressed  designs.  Six  normal  modes  of  vibration  were  used  in 
these  flutter  analyses  which  were  done  for  a Mach  number  of  1.6  and  an 
altitude  of  30,000  feet.  The  resulting  flutter  speeds  were  612  KEAS  and  624 
KEAS  for  the  unbalanced-  and  balanced-laminate  cases,  respectively. 
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Two  parallel  strength/flutter  redesign  studies  were  then  performed  to 
raise  the  flutter  speeds  of  both  designs  to  765  KEAS  (1.25  x 612).  Figure 
3.12  shows  the  history  of  the  resizing  steps.  Note  that  a flutter  spsed  close 
to  the  desired  value  was  attained  for  each  case  in  three  combined  resizing 
steps.  Two  additional  steps  were  performed  for  the  unbalanced  case  to 
approach  a converged  design. 

The  final  design  of  the  unbalanced-laminate  case  weighed  26  lb  more  than 
the  fully  stressed  design  from  which  it  was  derived.  Approximately  40%  of 
this  weight  increment  went  into  the  splice  plates  and  the  other  60%  was  distri- 
buted widely  among  the  composite  elements  in  various  regions  of  the  cover 
skins.  Layers  in  all  four  fiber  directions  were  affected.  Figure  3.13  shows 
the  unbalanced- laminate  upper  cover  layup  for  the  final  design.  A comparison 
of  this  layup  with  that  of  the  fully  stressed  design  (Figure  3.11)  reveals 
several  interesting  points,  which  are  discussed  below. 

First,  significant  amounts  of  material  were  added  to  the  metallic  splice 
plates  and  to  the  composite  material  in  their  vicinity.  For  the  composite 
material,  the  numbers  of  0°  fibers  aft  of  the  plates  and  the  ♦45°  and  90° 
fibers  forward  of  the  plates  (near  the  root  rib)  were  increased. 

Second,  material  was  added  in  the  middle  third  of  the  span  to  the  -45° 
fibers  and  even  more  so  to  the  ♦45°  fibers. 

Finally,  material  in  all  four  fiber  directions  was  added  to  the  Up  region, 
suggesting  the  presence  of  a mass-balance  effect  there,  rather  than  a stiff- 
ness effect.  Indeed,  inspection  of  the  fultter-velocity  derivatives  of  the 
variables  in  that  region  showed  that  the  derivatives  were  dominated  by  their 
kinetic-energy  components,  i.e.,  the  strain-energy  components  were  negligible. 
In  an  actual  design  effort,  it  would  be  more  practical  to  group  all  the  incre- 
mental tip  material  into  a single  mass-balance  variable  at  the  most  favorable 
location,  namely  at  the  tip  leading  edge.  Earlier  flutter  redesign  studies  on 
the  metallic  version  of  the  stabilizer  (see  Reference  1)  had  shown  that  mass 
balance  wac  extremely  effective  at  the  tip  leading  edge  - so  much  so  that  a 
near-optimum  design  had  been  obtained  with  a single  large  mass  at  the  . tip 
and  a relatively  small  amour',  of  cover  stiffening  near  the  root  rib  in  the 
vicinity  of  what  is  now  the  splice  plates.  This  approach  was  not  attempted  in 
the  current  study  because  a final  design  of  that  type  would  neither  test  nor 
demonstrate  the  major  new  capability  of  FASTOP-3.  amely  strength/flutter 
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r,,ilto.  composite  elements.  It  should  be  pointed  out,  hovever.  th.t  h.d 
sev.r.1  additionel  lultt.r  resizing  cycles  been  performed  in  the  current  stu  y. 
H Uh«.»  the,  the  program  «ou,d  h.ue  driven  the  design  to  one  in  uhtch 
mass-balance  was  the  dominant  contributor  to  increased  flutter  speed. 
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